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INTRODUCTION 


Breast  cancer  has  been  a  major  cause  for  death  for  women  in  the  United  States. 
Currently,  the  major  tool  for  screening  and  diagnosing  breast  cancer  is  the  X-ray 
mammography.  With  this  technique.  X-ray  images  of  breasts  at  low  kVp  to  produce 
reasonably  high  image  contrast  for  detecting  microcalcifications,  masses  and  examining 
tissue  structure  of  the  breast.  Because  microcalcifications  are  generally  small  in  size, 
high  spatial  resolution  is  required  to  image  them  with  reasonable  sharpness  and  contrast 
for  detection.  On  the  other  hand,  high  detection  efficiency  is  required  to  minimize  X-ray 
exposure  and  therefore  cancer  risk  to  the  breast  while  producing  reasonable  signal-to- 
noise  ratio  for  detecting  the  subtle  contrast  of  small  calcifications  as  well  as  tumor 
masses.  Current  mammographic  technique  employs  high  resolution  screen-film 
combination  as  the  image  detector  and  recording  system.  This  system  provides  adequate 
spatial  resolution  and  detection  efficiency  for  the  mammographic  imaging  tasks. 

Although  screen-film  techniques  have  provided  reasonable  image  quality  for 
various  applications,  there  are  several  drawbacks.  First  of  all,  screen-film  combinations 
require  accurate  control  of  X-ray  exposures  to  avoid  over-  or  under-exposure.  Thus, 
retakes  due  to  incorrect  exposures  represent  a  significant  ratio  of  all  radiographs  taken. 
This  increases  film  consumption  as  well  as  patient  exposures.  Secondly,  film  images  are 
typically  displayed  on  viewing  boxes  with  fixed  contrast  and  brightness.  They  can  not  be 
enhanced  or  processed  without  first  being  digitized.  Lastly,  storage,  retrieval  and 
transport  of  films  are  becoming  an  increasingly  serious  problem  when  films  have 
accumulated  for  years. 

With  recent  advances  of  digital  technology,  the  Picture  Archival  and 
Communiction  Systems  have  become  technically  and  economically  feasible  for  use  in  the 
radiology  environment.  Such  systems  are  expected  to  significantly  improve  the 
management  of  radiological  images.  However,  with  such  systems  in  use,  it  also  becomes 
more  important  to  have  a  high  quality  but  economic  technique  for  direct  digital 
acquisition  of  projection  radiography  images.  Various  digital  radiography  techniques 
have  been  developed  to  acquire  projection  radiography  images  .  Some  of  them  went 
further  to  provide  improved  image  quality  by  incorporating  scatter  rejection  capability  or 
using  X-ray  detectors  of  superior  quantum  efficiency.  However,  many  of  these 
techniques  also  are  either  too  expensive  to  implement  or  suffer  fiom  poor  spatial 
resolution.  Among  these  digital  radiography  techniques,  the  storage  phosphor  imaging 
technique  has  evolved  to  be  a  promising  candidate  for  large  scale  implementation  of 
digital  radiography  capability  in  radiology  enveironment. 

Storage  phosphor  (SP)  imaging  is  a  technique  of  using  metastable  electron  states 
to  temporarily  store  X-ray  exposure  information  and  then  using  scaiming  laser  beam  to 
read  out  the  latent  image  signal  as  light  emitted  when  the  electrons  in  the  metastable  state 
are  excited  and  return  to  the  ground  state.  Because  the  SP  image  plates  can  be  housed 
and  exposed  in  cassettes  similar  to  those  used  in  conventional  screen-film  radiography. 
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the  SP  imaging  technique  can  be  implemented  without  replacing  the  currently  used  X-ray 
equipment  and  with  minimal  training  requirement  for  the  technologists.  Thus,  it  provides 
a  convenient  and  economic  way  to  incorporate  digital  radiography  capability. 

Although  the  SP  imaging  technique  has  been  actively  used  in  foreign  countries  for 

mammographic  imaging  applications,  its  applications  in  the  United  States  have  been 

limited  to  general  radiography,  particularly  to  chest  imaging  in  the  intensive  care  unit 

Due  to  technical  or  marketing  considerations,  typical  commercial  SP  imaging  systems 

have  been  configured  to  use  a  pixel  size  ranging  from  86  to  200  pm  with  a  laser  spot  size 

8-10 

of  100  pm  or  greater.  Due  to  this  limitation  of  the  screen  optics  ,  current  commercial 
CR  systems  do  not  have  sufficient  spatial  resolution  capability  for  mammographic 
imaging  applications  ’  .  submits  the  spatial  resolution  capability  for  CR  imaging.  This 
practice  poses  limitations  on  the  use  of  current  commercial  SP  imaging  systems  for  high 
resolution  imaging  applications,  among  which  is  the  mammographic  imaging. 


Overall  Goals 


In  this  project,  we  propose  to  modify  and  optimize  an  experimental  CR  system  for 
high  resolution  mammographic  imaging  applications.  Image  properties  of  the  optimized 
system  will  be  characterized  by  measuring  the  signal  response,  MTF  and  noise  properties. 
A  eomprehensive  comparison  study  based  on  phantoms  will  be  conducted  to  evaluate  and 
compare  the  CR  technique  with  the  conventional  screen-film  technique  for  detection  of 
microcalcifications  as  well  as  low  contrast  soft  tissue  masses.  Possibility  for  patient  dose 
reduction  will  also  be  investigated.  In  an  addendum  to  our  proposal  dated  October  26, 
1992,  we  also  proposed  to  perform  preliminary  investigations  on  data  compression  of 
mammographic  or  similar  types  of  images. 

Main  Tasks 


The  project  consists  of  five  major  components:  (1)  system  modification  and 
optimization,  (2)  system  characterization  and  quality  control,  (2)  comparison  study  based 
on  phantom  images,  (4)  comparison  study  based  on  ROC  studies  and  (5)  image 
compression  study.  These  tasks  have  largely  been  accomplished  and  will  be  described 
and  discussed  in  details  in  the  following  sections. 
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1.  SYSTEM  MODIFICATION  AND  OPTIMIZATION 


Theoretical  Considerations 

In  storage  phosphor  imaging,  the  X-rays  expose  the  SP  image  plates,  leaving  a 
two-dimensional  distribution  of  trapped  electron  energy.  This  energy  distribution  is  also 
referred  to  as  the  latent  image  signal.  During  image  readout,  the  trapped  electron  energy 
is  released  by  a  scanning  laser  beam  ’  ’  .  The  blurring  processes  for  the  SP  image  readout 
is  illustrated  in  Figure  1.  The  latent  image  signals  is  generated  within  a  few  microns  from 
where  the  X-ray  photon  interact  with  the  screen  material.  Therefore,  there  is  minimal 
loss  of  spatial  resolution  in  this  process.  However,  the  laser  beam  has  a  finite  size  which 
may  be  characterized  by  an  intensity  distribution  at  entrance  surface  of  the  screen.  In 
addition,  the  laser  beam  is  subject  to  photon  scattering  and  light  diffusion  once  it  enters 
the  screen^’^.  Both  facts  result  in  in  a  degradation  of  the  spatial  resolution  quality*’". 
Although  the  propagation  of  the  laser  beam  inside  the  screen,  the  excitation  of  the  trapped 
electrons  and  release  of  their  energy  are  rather  complicated  phenomena,  the  blurring 
process  can  be  described  by  a  simple  model.  In  this  model,  the  laser  beam  profile  and 
light  diffusion  in  the  screen  are  characterized  by  the  point  spread  fimctions  (PSF’s),  PSFg 
and  PSFs.  The  blurring  process  may  be  expressed  as  follows: 

S{x,  y)  =  k^^dx''  dy'b  (x" ,  y" )  j|  dx'  dy'  PSF^  (x"-x' ,  y'-y'  )lg(x'-x,  y'-y) 

=  kP,  jldx"dy'b(x",y')l\dx'dyPSF,ix"-x\y"-y')PSF,ix'-x,y-y) 


where 


PSFs(x,y) 


^^dx' dy' lyp)  Pq 


Ig{x'-x,  y'-y)  =  laser  beam  intensity  profile  at  (x, 3^)  in  the  detector  plane, 
PSFg  (x"-x' ,  y"-y )  =  point  spread  function  for  laser  light  diffusion  at  (x' ,  y' ), 
a  (x"  ,y")=  area  density  of  trapped  electron  energy. 

a(x",>'")  and  S{x,y)  may  be  viewed  as  input  and  output  images  respectively. 


The  above  equations  may  be  alternatively  rewritten  as  follows 


S  =  kP^-  PSFg  0  (PSFs  ®  cr ) 

=  kPo-PSF0a 

and  the  pre-digitization  point  spread  function,  PSF,  is  formed  from  PSF^  and  PSFs  as 
follows: 
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X-rays 


Scanning 


Figure  1  Schematic  drawing  for  the  image  blurring  processes  in  CR  imaging 
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PSF  =  PSFj,0PSFs 

PSF,  PSFg  and  PSFg  may  be  used  to  compute  the  line  spread  functions  (LSF’s)  and 
MTF’s.  The  resulting  pre-digitization  MTF,  MTF,  is  related  to  the  beam  profile  and 
screen  MTF’s,  MTFg  and  MTFs,  as  follows: 

MTF  =  MTFg  0  MTFg 

From  the  above  equation,  it  is  clear  that  both  the  beam  profile  and  screen  MTF’s  are 
important  in  determining  the  overall  MTF  or  the  spatial  resolution  quality  of  the  imaging 
system.  A  large  beam  size  or  poor  screen  MTF  or  both  would  all  degrade  the  resolution 
quality  of  the  imaging  system. 

To  estimate  the  effect  of  the  beam  spot  size  on  the  MTF,  we  may  assume  that  the 
beam  profile,  is  a  Gaussian  function,  with  the  full  half  maximum  width,  .  Based 

on  this  assumption,  the  point  spread  function,  line  spread  function  and  MTF  can  be 
derived  to  be  as  follows: 

.2 

-41n2~ 

41n2 

- T'^  ' 

ndi 


LSF,{x,y)=  \dyPSF,{x,y) 

x2 

-4  In  2  — j 

2Vin2 

2 

— - - 

MTF,if)  =  e  41n2  ^ 

Using  the  above  equation,  the  beam  profile  MTF  is  computed  and  plotted  as  a  function  of 
the  spatial  frequency  for  various  beam  sizes  {d^  =  22,  43,  65  and  86  pm)  in  Figure  2.  At 

2,  5  and  10  cycles/mm,  the  beam  profile  MTF  values  for  a  beam  size  of  43pm  are  97.4, 
84.8  and  51.8%  respectively.  This  means  that  with  a  beam  spot  size  of  43  pm,  only 
slightly  over  50%  of  the  screen  MTF  at  10  cycles/mm  is  preserved  in  the  overall  MTF 
while  those  at  2  and  5  cycles  are  mostly  preserved.  A  switch  to  22  pm  spot  size  would 
increase  the  MTF  at  10  cycles/mm  to  84.8%  and  preserve  most  of  the  screen  MTF  in  the 
overall  MTF  at  frequencies  of  up  to  10  cycles/mm.  On  the  other  hand,  an  increase  of  the 
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Frequency  (cycles/mm) 


Figure  2  MTF’s  computed  for  four  different  laser  beam  spot  sizes  (half  maximum 
width  =  0.022,  0.043,  0.065  and  0.086  mm)  for  frequencies  of  up  to  12 
cycles/mm.  The  beam  profiles  are  assumed  to  have  the  shape  of  Gaussian 
functions. 
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beam  spot  size  to  86p,m  would  decrease  the  beam  profile  MTF  to  7.2%.  Combined  with 
an  already  low  screen  MTF,  this  would  reduce  the  overall  MTF  to  almost  zero  level. 

System  Modifications 

Figure  3  shows  a  schematic  diagram  of  the  laser  scanning  optics  used  in  our 
experimental  image  reader.  The  collimated  beam  from  an  argon  laser  unit  passes  through 
a  lens  A  and  focuses  at  the  acoustic  optical  modulator  (AOM).  The  divergent  output  is 
collimated  again  by  lens  C,  reflected  by  mirror  D  and  focused  again  by  a  microscopic 
objective  lens  F.  An  aperture  E  is  used  to  block  the  zero  order  peak  and  select  the  first 
order  peak  of  the  AOM  output.  The  latter  has  an  intensity  adjustable  through  digital 
control  of  the  AOM.  The  output  of  the  lens  F  is  collimated  into  a  widened  beam  by  lens 
G  (beam  expander)  which  is  reflected  by  the  mirror  H  and  deflected  by  the  digitally 
controlled  galvanometer  mirror  I.  The  deflected  beam  passes  through  an  F-theta  or  laser 
scan  lens  J  and  focuses  on  the  image  plane  K. 

The  laser  spot  size  at  the  image  plane  K,  d,  is  determined  by  many  factors 
including  the  focal  lengths  of  the  lenses  A,  E  and  J.  The  diameter  of  the  expanded  beam 
at  the  output  of  the  lens  G  and  the  size  of  the  galvanometer  mirror  I  also  pose  limit  on 
how  small  the  laser  spot  size  can  be.  In  1993,  we  successfully  reduced  the  laser  spot  size 
from  120fim  to  SOpm  by  changing  the  focal  length  of  lens  A  from  250mm  to  125nim  and 
replacing  Lens  F  (previously  a  5x  microscopic  objective  lens)  with  a  lOx  one.  As 

discussed  before,  the  l/e^  diameter,  c,  is  often  used  to  describe  the  beam  size  in  laser 
optics.  (7  can  be  estimated  from  d  as  follows: 


Thus,  the  l/e^  diameter  was  reduced  from  204pm  to  approximately  85pm. 

To  further  reduce  the  beam  spot  size,  we  also  replaced  the  AOM( Acoustic-Optical 
Modulator)  with  a  more  efficient  one  and  move  it  to  between  the  laser  unit  and  Lens  A. 
This  is  possible  because  the  new  AOM  does  not  required  focused  beam  and  can  operate 
on  direct  output  from  the  laser  imit.  A  small  aperture  (10-3  0pm  in  diameter)  will  be 
positioned  at  the  focal  point  between  Lens  A  and  Lens  C.  The  use  of  this  aperture  will 
force  the  beam  spot  size  to  be  reduced  throughout  the  rest  of  the  optics.  However,  a 
significant  amount  of  laser  power  will  be  sacrificed  through  blocking  the  light  outside  the 
aperture. 

Another  source  for  loss  of  laser  power  is  through  the  beam  expander  G, 
galvanometer  mirror  I  and  laser  scan  lens  J.  During  our  first  stage  of  system 
modification,  adequate  but  marginal  laser  power  was  achieved  with  currently  used  beam 
expander,  galvanometer  mirror  and  laser  scan  lens.  However,  they  could  limit  the  laser 
power  to  unusable  range  if  further  reduction  of  the  beam  spot  size  (down  to  25-3 0pm)  is 
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I  :  ga I vanometer 


Figure  3  Schematic  drawing  for  the  laser  scanning  optics  in  an  experimental  CR  image 
reader. 
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desired.  Instead  of  replacing  these  components,  a  easier  approach  is  to  use  a  more 
efficient  AOM  (Multiwavelength  Visible  Light  Modulator,  Model  No.  N48062-2.5-.55, 
NEOS  Technologies  Inc.)  and  increase  the  power  available  at  the  beginning.  The  new 
AOM  has  a  maximum  light  transmission  of  over  90%  and  therefore  can  produce  a  much 
higher  power  than  the  previously  used  one  (made  by  Newport  Research  Inc.,  with  an 
efficiency  of  less  than  10%).  This  higher  power  provides  us  a  greater  flexibility  in  further 
reducing  the  beam  spot  size. 

We  also  acquired  a  number  of  high  resolution  storage  phosphor  screens  for  use  in 
image  acquisition  for  the  ROC  studies.  These  screens  include  the  newest  one  from  Fuji: 
HR-V.  However,  our  test  showed  that  HR-V  offers  a  spatial  resolution  quality  similar  to 
those  of  the  HRIII  series  although  it  does  produce  a  higher  light  gain. 


Measurement  of  laser  beam  soot  size 


As  discussed  in  the  Theory  section,  laser  beam  spot  size  is  one  of  the  major 
factors  that  dictate  the  spatial  resolution  quality  of  the  resulting  CR  images.  To  evaluate 
the  laser  spot  size,  two  optical  slit-photo  Detector  assemblies  were  constructed  and 
attached  to  one  end  of  the  image  plate  carrier  stage.  Each  assembly  consists  of  a  5  pm 
optical  slit  and  a  2.5mm  slit  mounted  in  front  of  a  photo-transistor  (Figure  4).  The  photo 
transistor  output  is  amplified  and  then  forwarded  to  an  oscilloscope  for  monitoring.  It 
can  also  be  digitized  for  analysis  in  the  computer.  However,  a  dedicated  high  speed 
analog-to-digital  converter  would  be  required  to  measure  sub-millimeter  spot  sizes  with 
reasonable  precision.  This  option  is  being  planned  but  not  implemented  yet. 

To  measure  the  laser  beam  spot  size,  the  assemblies  were  first  positioned  in  the 
laser  scanning  path.  As  the  laser  beam  scarmed  across  the  two  optical  slits  in  each 
assembly,  a  signal  consisting  of  two  pulses  was  generated  (Figure  5).  This  signal  is 
proportional  to  the  laser  beam  power  profile,  integrated  in  the  longitudinal  direction  (in 
which  the  translational  stage  moves).  The  slit  signal  can  be  observed  on  a  scope  to 
measure  the  full  width  at  half-maximum  of  the  two  pulses.  Assuming  the  (half 
maximum)  pulse  widths  for  the  5pm  and  2.5mm  slits  are  t  and  T  respectively,  the  half 
maximum  width  of  the  laser  beam  profile  on  the  image  plane,  d,  can  be  estimated  as 
follows: 


d  =  2.5mm  x 


The  process  of  measuring  the  laser  beam  spot  size  was  later  automated.  Instead  of 
being  displayed  on  a  scope,  the  slit  signal  is  digitized  and  analyzed  in  a  microcomputer  to 
provide  both  the  half-maximum  and  widths.  The  measurement  can  be  repeated  at  a 
rate  of  approximately  1  per  second  to  allow  the  optics  to  be  adjusted  while  monitoring  the 
beam  spot  size.  The  difficulty  with  measuring  the  laser  spot  size  resides  with  the  location 
of  the  maximum  and  background  levels.  Errors  in  determining  these  levels  may  bias  the 
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measurements.  With  continuing  refinements,  the  smallest  beam  spot  size  we  could  have 
achieved  was  measured  to  be  30  pm  for  the  FWHM  or  of  70pm  for  e'  width  in  the  high 
resolution  mode.  This  corresponds  to  a  MTF  value  of  approximately  70%  at  10 
cycles/mm. 

Further  reduction  is  possible.  However,  it  may  start  to  degrades  the  sampling 
process  and  generate  gaps  where  x-ray  information  would  be  lost  and  image  noise  may 
increase.  Although  we  continue  to  investigate  the  pros  and  cons  of  smaller  spot  size,  we 
have  reached  the  conclusion  that  based  on  practical  considerations,  it  would  be  a  good 
compromise  to  use  40pm  laser  beam  spot  in  our  comparison  study. 

We  have  also  developed  a  series  of  techniques  for  testing  the  operation  of  ESPIS. 
One  of  these  techniques  is  direct  measurement  of  the  scanning  laser  intensity  pattern.  We 
have  been  able  to  use  film  to  monitor  and  measure  the  spatial  intensity  distribution  of  the 
scanning  laser  beam.  This  technique  has  proved  to  be  useful  in  checking  the  stability  of 
the  laser  power  or  the  proper  function  of  the  AOM. 

System  Configuration 

The  CR  imaging  system  used  in  this  study  is  an  experimental  system  built  on  an 
optical  bench  with  standard  optics  components  and  a  high  power  argon  laser  unit.  The 
system  employs  galvano-mirror  to  provide  a  scanning  laser  beam.  During  scan  readout, 
the  SP  image  plate  is  placed  on  a  translational  stage  and  moved  into  the  scanning  laser 
beam  for  signal  readout.  The  system  is  capable  of  scanning  and  digitizing  with  a  variety 
of  pixel  sizes,  including  the  frequently  used  172,  86  and  43  pm.  The  readout  signals  are 
digitized  into  12  bit  data  without  going  through  logarithmic  amplification  or  mapping. 
An  adjustable  active  bandpass  filter  system  is  used  to  match  the  bandpass  frquency  of  the 
signal  input  to  the  analog-to-digital  converter. 

A  LoRAD  M-III  mammogrpahic  unit  was  used  to  provide  X-ray  exposures.  A 
variety  of  image  plates  are  available  for  use  with  our  SP  system.  For  this  study,  Fuji  HR 
Illn  and  ST  Illn  plates  (cut  to  size  of  mammographic  screens)  were  housed  in  Kodak 
Min-R  cassettes  for  exposures.  This  unit  employs  molybdenum  target  and  reciprocating 
grid.  For  screen-film  images,  Kodak  Min-R  screens  are  used  in  conjimction  with  Kodak 
Min-R  E  film. 
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SYSTEM  CHARACTERIZATION 


Modulation  Transfer  Function  (MTF) 

We  continued  to  use  a  lead  bar  resolution  pattern  (Nuclear  Associates  07-553)  to 
monitor  and  measure  the  spatial  resolution  quality  of  the  ESPIS.  It  consists  of  22  groups 
of  bar  patterns  with  a  frequency  ranging  from  0.25  to  10  Ips/mm.  Each  group  has  a  4 
cycle  bar  pattern.  A  specification  list  for  the  resolution  bar  pattern  is  attached  as 
Appendix  I. 

The  bar  pattern  images  are  visually  examined  to  evaluate  the  spatial  resolution 
limit  of  the  imaging  system.  A  fully  resolved  pattern  (with  good  contrast)  generally 
indicates  that  the  MTF  is  around  0.05  or  greater  at  that  particular  frequency.  A  partially 
or  barely  resolved  pattern  with  faint  contrast  indicates  that  the  MTF  is  probably  between 
1%  and  5%.  Unresolved  bar  patterns  generally  indicate  that  the  MTF  is  below  1%.  This 
observation  could  vary  from  person  to  person  but  may  be  used  a  crude  guideline  when 
visually  examining  the  bar  pattern  images. 

The  methods  for  measuring  the  spatial  resolution  properties  of  a  storage  phosphor 
imaging  system  have  been  previously  discussed  and  demonstrated^’9’12,13  Our 
methods  focus  on  the  use  a  standard  resolution  test  pattern  as  described  before.  The 
square  wave  response  function  (SWRF)  can  be  measured  by  plotting  the  signal  profile 
across  the  bar  patterns  and  then  measuring  the  average  amplitude  for  each  frequency. 
Because  the  square  waves  contain  high  order  harmonics,  the  SWRF  can  not  be  directly 
used  to  compute  the  MTF  without  correction  for  these  harmonics.  We  have  devised  and 
tested  a  signal  processing  method  to  compute  the  MTF  from  the  signal  profile  of  a  bar 
pattern. 


The  MTF  is  measured  from  the  digital  SP  images  of  a  resolution  bar  pattern 
(Nuclear  Associates  07-553)  to  evaluate  the  spatial  resolution  quality  of  the  resulting 
images.  This  pattern  provides  bar  patterns  at  frequencies  ranging  from  0.25  to  10 
cycles/mm.  To  measure  the  MTF,  the  image  of  the  bar  patterns  is  first  segmented  into 
sections  of  various  frequencies.  Signal  profiles  across  the  bar  patterns  from  each  section 
are  then  concantenated  over  20  lines  into  a  long  array  of  signals  modulated  at  a  single 
frequency.  This  array  is  offset  to  correct  for  x-ray  transmission  through  lead  bars  and  the 
presence  of  scatter  signals.  The  Fourier  transform  of  this  long  array  is  then  computed  to 
and  used  to  obtain  a  frequency  spectrum  of  the  signals.  This  spectrum  contains  a  primary 
peak,  corresponding  to  the  frequency  of  bar  patterns,  and  higher  order  harmonics  at 
multiples  of  this  frequency.  These  harmonics  are  ignored  in  the  computation  of  the  MTF. 
The  height  of  the  primary  peak,  S(f),  is  normalized  by  the  DC  component,  S(0),  to 
compute  the  MTF  at  the  frequency  at  f  as  follows: 
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MTF^ 


S{f) 

0.6369-5(0) 


where  the  correction  factor  0.6369  is  the  height  of  the  primary  peak,  normalized  by  the 
DC  component,  in  a  perfect  square  wave,  toggling  between  0  and  some  positive  value. 
This  method  provides  an  rather  easy  way  to  evaluate  the  MTF  or  spatial  resolution 
properties  of  a  digital  imaging  system. 

We  have  measured  and  compared  the  system  MTF’s  for  two  different  laser  beam 
spot  sizes:  30  and  50pm  (FWHM).  The  results,  shown  in  Figure  6,  indicate  a  significant 
improvement  of  the  MTF.  Notice  that  the  ratio  of  improvement  is  rather  insignificant  at 
low  frequencies.  However,  it  becomes  more  and  more  significant  as  the  frequency 
increases.  The  MTF  value  almost  doubles  at  frequencies  above  6  cycles/mm.  Despite  the 
improvement,  the  MTF  values  are  still  low  at  high  frequencies.  Further  increase  of  the 
MTF  would  require  a  significant  improvement  of  the  CR  screens. 

We  have  also  measured  the  MTF  as  a  function  of  the  laser  beam  power,  ranging 
from  8  to  12  mW.  The  results,  shown  in  Figure  7,  indicate  that  the  MTF  does  not  depend 
on  the  laser  beam  power  in  the  useful  range  of  8-12  mW.  It  should  be  noted  that  outside 
this  range,  the  power  may  be  too  strong  to  avoid  overheating  or  even  burning  the 
phosphor  or  too  low  for  the  position  encoding  to  properly  function. 

In  Figure  8,  the  MTF  is  plotted  as  a  function  of  the  frequency  for  both  the  HR  and 
ST  screens  to  demonstrate  the  improvement  from  using  a  HR  screen.  The  MTF  with  the 
HR  screen  is  generally  higher  at  all  frequencies  except  at  very  low  (near  0.25  cycles/mm) 
and  very  high  frequency  (near  10  cycles/mm).  The  closeness  of  the  HR  and  ST  MTF 
values  near  10  cycles/mm  is  rather  puzzling.  However,  it  may  be  partially  explained  by 
the  fact  that  in  mammographic  imaging  more  X-rays  are  absorbed  near  the  screen  surface 
(assuming  that  X-rays  enter  the  screen  from  the  surface)  and  high  frequency  signals  may 
be  mostly  generated  near  the  screen  surface. 

In  Figure  9,  the  MTF  of  our  high  resolution  CR  system  is  compared  with  a  typical 
commercial  system  operated  in  two  different  techniques  (an  early  KESPR  prototype  by 
Eastman  Kodak  Company).  One  technique  employs  172pm  pixel  size  to  generate  2048  x 
2500  images  for  14”xl7”  field  size.  This  is  configured  for  medium  resolution 
applications  like  chest  imaging.  A  second  typical  technique  employs  100pm  pixel  size  to 
generate  2048  x  2500  images  for  the  smaller  8”xl0”  field  of  view.  This  configuration 
provides  slightly  better  resolution  and  is  intended  for  more  demanding  applications  like 
pediatric  or  orthopedic  imaging.  Figure  9  shows  that  our  high  resolution  CR  system 
produces  significantly  better  MTF  than  either  techniques  with  the  commercial  system. 
Because  similar  CR  screens  were  used  (ST  for  14”xl7”  and  HR  for  8”xl0”),  the  lower 
MTF’s  of  the  two  commercial  techniques  are  probably  due  to  the  larger  spot  size  used. 
The  larger  pixel  sizes  used  by  the  commercial  systems  may  also  contribute  to  the 
degradation  of  the  MTF’s. 
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Figure  6  MTF’s  of  a  CR  imaging  system:  30|am  versus  50|j.m  laser  beam  spot  sizes. 
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MTF  of  FJ  plates,  LP:  8(*).  10(+)  12(x)and  16(o)  images 


Figure  7  Effect  of  laser  power  on  the  MTF.  The  MTF  was  measured  for  various  laser 
powers,  ranging  from  8  to  16  mW:  8  mW,  10  mW,  "x";  12  mW  and 

"o":16mW. 
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Figure  8 


Frequency  (cycles/mm) 


MTF  of  the  high  resolution  CR  imaging  system:  HR  screen  versus  ST 
screen. 
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In  Figure  10,  the  MTF  of  our  high  resolution  CR  system  measured  with  a 
mammographic  exposure  (25  kVp,  10  mAs)  is  compared  with  that  measured  with  a 
typical  chest  exposure  (110  kVp,  6  mAs).  The  plots  show  that  the  MTF  is  significantly 
higher  with  mammographic  exposures.  This  may  be  explained  by  the  fact  that  high 
frequency  signals  are  mostly  generated  near  the  surface  of  the  screen  and  the  latent  image 
signals  are  distributed  more  toward  the  screen  surface  with  the  msimmographic  X-rays 
than  with  chest  X-rays.  MTF’s  measured  with  high  kVp  mammographic  exposures 
should  be  between  the  two  MTF’s  shown  in  Figure  10. 

Films  1  shows  a  screen-film  image  of  a  resolution  test  pattern  (Nuclear 
Associates,  Model  07-553)  acquired  with  a  mammographic  unit  at  25  kVp  and  6  mAs. 
Films  2,  3  and  4  show  CR  images  of  the  same  pattern  acquired  with  three  different 
techniques:  (Film  2)  exposed  with  a  ST  screen  and  read  out  with  30pm  laser  beam  size, 
(Film  3)  exposed  with  a  HR  screen  and  read  out  with  30pm  laser  beam  size  and  (Film  4) 
exposed  with  a  HR  screen  and  read  out  with  50pm  laser  beam  size.  All  images  were  read 
out  with  a  pixel  size  of  43  pm.  Notice  that  the  screen-film  image  displays  best  resolution: 
all  bar  patterns  are  clearly  resolved.  The  image  in  Film  3,  acquired  with  HR  screen  and 
small  beam  spot  size,  shows  the  best  spatial  resolution  quality  among  the  three  CR 
images.  Large  beam  spot  size,  thicker  ST  screen  or  both  degrade  the  spatial  resolution 
quality  of  the  imaging  system. 

Detective  Quantum  Efficiency 

Detective  quantum  efficiency  (DQE)  is  often  measured  as  a  function  of  the  spatial 
frequency  and  used  to  characterize  the  capability  of  an  imaging  system  to  preserve  the 
signal-to-noise  ratio  (SNR)  in  the  input  image  signals.  DQE  is  generally  defined  as 
follows: 


DQE  = 


SNRJ 


where  S7Vi?,„  and  are  the  image  signal-to-noise  ratios  at  the  input  and  output.  To 

measure  the  DQE  as  a  function  of  the  spatial  frequency,  f  (Ips/mm),  we  used  the 
following  expression  for  the  frequency  dependent  DQE: 


DQEif)  = 


SNR.. 


[nps(J)) 


■MTF{ff 


where  NPS(f)  is  the  noise  power  spectrum,  S  is  the  average  signal  over  the  area  where 
NPS(f)  is  measured,  MTF(f)  is  the  modulation  transfer  function.  This  equation  is  similar 
to  those  previously  used  in  literature  on  DQE  measurement^’9,1 1,15 
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Figure  10  MTF’s  of  the  high  resolution  SP  imaging  system:  mammographie  technique 
versus  chest  imaging  technique. 
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We  have  measured  the  noise  power  spectrum  (NPS)  and  detective  quantum 
efficiency  (DQE)  for  our  CR  system.  To  compute  the  DQE(t),  NPS(f)  is  computed  from 
image  data  over  an  area  with  largely  uniform  signals  while  SNR^^  is  computed  from  the 
number  of  photons  which  is  estimated  from  the  x-ray  techniques  used.  Figure  1 1  shows 
the  noise  power  spectrum  computed  from  the  image  data  in  a  largely  uniform  region.  The 
x-ray  techniques  used  in  acquiring  the  image  data  are  25  kVp,  10  mAs,  with  grids  at  an 
SID  of  65  cm.  The  resulting  detector  exposure  is  approximately  6  mR.  Notice  that 
because  of  image  blurring,  the  noise  spectrum  is  no  longer  white.  Instead,  it  decreases 
with  increasing  frequency.  Figure  12  shows  the  DQE  plotted  as  a  function  of  the 
frequency. 

The  reduction  of  laser  beam  spot  size  affects  the  DQE  mainly  through  the  MTF 
improvement.  However,  if  the  laser  beam  size  is  too  small  as  compared  to  the  pixel  size, 
loss  of  quantum  detection  efficiency  may  occur  and  degrade  the  DQE  at  all  frequencies. 


Quality  Control 

Quality  control  of  the  CR  imaging  system  has  been  performed  by  acquiring  and 
visually  checking  images  of  the  resolution  bar  pattern,  CIRS  phantom,  RMI  phantom  and 
the  specially  designed  aluminum  wire  phantom.  Because  some  improvements  have  been 
made  on  the  system,  the  criteria  for  visual  check  has  slightly  changed  now.  In  general, 
for  the  resolution  bar  pattern  image,  we  expect  to  see  and  resolve  all  bars  at  frequencies 
of  up  to  8.5  Ips/mm  fully  and  clearly.  The  10  Ips/mm  bars  should  also  be  fully  resolved 
but  with  subdued  contrast.  Partially  resolution  of  the  10  Ips/mm  bars  is  an  indication  of 
presence  of  structural  artifacts.  The  images  of  the  CIRS  phantom  provides  a  convenient 
semi-quantitative  method  for  checking  the  overall  quality  of  the  system,  including  the 
spatial  resolution  and  noise  properties.  Same  technique  (25kVp,  144mAs,  65cm  SID, 
with  antiscatter  grids)  will  be  used  in  all  quality  control  exposures.  In  general,  we  expect 
to  see  all  microcalcifications  with  a  diameter  of  200|j,m  or  greater  and  all  masses  with  a 
diameter  of  2nim  or  greater.  All  fibrils  should  be  fully  resolved.  The  RMI  phantom 
offers  a  more  qualitative  check.  All  images  should  have  similar  apparent  quality  as 
previous  images.  The  aluminum  wire  phantom  provides  a  finer  technique  to  grade  the 
system  performance.  At  25  kVp  and  26  mAs,  we  expect  to  resolve  the  1 "  long  wires  as 
thin  as  50)am  in  diameter  and  small  wire  segments  as  thin  as  125pm. 

Quality  control  of  mammographic  X-ray  units  and  the  screen-film  combinations 
have  been  previously  discussed  ^  1 04 1 1 .  The  quality  of  the  mammographic  X-ray  unit  is 
routinely  monitored  and  maintained  by  the  clinieal  staff 
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Frequency  (cycles/mm) 


Figure  1 1  Noise  power  spectrum  (NPS)  computed  from  image  data  in  an  area  of  largely 
uniform  brightness.  The  spectrum  was  normalized  by  dividing  it  by  the  square 
of  the  average  signal. 
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COMPARISON-  PHANTOM  IMAGES 


Phantom  Selection  and  Design 

Image  quality  of  a  mammographic  unit  is  often  checked  with  phantoms  simulating 
a  breast  with  the  three  major  disease  symptoms:  microcalcifications  and  masses.  Various 
phantoms  have  been  designed  for  this  purpose.  We  have  used  two  commercial  phantoms 
to  test  the  image  quality  and  compare  different  imaging  systems  or  technical  factors.  The 
first  one  is  a  Tissue-Equivalent  Breast  Phantom  made  by  the  Computerized  Imaging 
Reference  Systems  (CIRS),  Inc.  of  Norfolk,  Virginia.  It  consists  of  a  breast  shape 
phantom  made  of  tissue  equivalent  material,  with  a  step  wedge,  groups  of 
microcalcifications,  masses  and  fibrils  embedded  inside  the  phantom.  A  specification 
sheet  for  the  phantom  is  attached  as  Appendix  II.  The  step  wedge  consists  a  square 
region  with  a  contrast  corresponding  to  additional  5mm  of  glandular  tissue  over  the 
surrounding  background  tissue  and  a  second  square  region  with  a  contrast  corresponding 
to  additional  5mm  of  fat  tissue.  The  12  microcalcification  groups  ranges  from  0.12mm  to 
0.35mm  in  diameter.  The  6  masses  range  from  1mm  to  6mm  in  diameter.  The  fibrils 
have  a  diameter  of  8.7p,m.  This  phantom  provides  a  good  tool  for  comparing  the  overall 
image  quality  of  different  imaging  systems,  techniques  or  technical  factors. 

The  second  phantom  is  the  RMI  Anthropomorphic  Breast  Phantom  (Model  165, 
Gammex-RMI).  The  phantom  consists  of  a  breast  tissue  equivalent  (50%  adipose,  50% 
glandular)  plastic  casting  and  a  high  resolution  detail  layer  enclosed  in  a  protective 
acrylic  case.  Together,  these  components  produce  radiographic  images  that  are  similar  to 
a  true  mammogram..  To  aid  the  evaluation  of  system  performance,  the  phantom  also 
produces  a  nine  level  gray  scale  step  wedge  and  a  group  of  resolution  line  pairs  in  the 
mammogram.  The  phantom  is  unique  in  that  the  characteristics  of  its  components  were 
derived  from  an  actual  mammogram.  The  three  dimensional  surface  structure  of  the 
breast  tissue  equivalent  plastic  casting  produces  the  low  and  medium  resolution  image 
detail.  High  resolution  image  details  are  produced  by  an  attenuating  coating  of  a  stable 
mercury-silver  amalgam  on  photographic  film.  In  addition  to  protecting  the  contents,  the 
acrylic  case  acts  as  a  uniform  attenuation  layer  across  the  entire  phantom 

The  advantage  of  using  the  above  two  phantoms  is  that  they  more  or  less  resemble 
a  real  breast  in  tissue  thickness,  composition  and  even  structures.  Thus  they  allow  the 
performance  to  be  checked  in  a  realistic  fashion.  However,  they  are  more  suitable  for 
quality  control  than  for  critical  comparison  of  image  quality.  The  simulated 
microcalcifications  in  the  CIRS  phantom  have  a  size  of  at  least  120p.m.  In  addition,  the 
sizes  of  different  groups  increment  by  40pm.  To  facilitate  more  critical  and  quantitative 
comparison  of  image  quality,  we  have  developed  a  phantom  referred  to  as  the  perception 
phantom.  The  perception  phantom  was  constructed  of  aluminum  wires  of  8  different 
sizes:  2,  3,  4,  5,  6,  7,  8  and  16  mils.  These  wires  were  cut  into  about  1"  long  and  placed 
in  radial  directions  over  a  1 "  thick  Lucite  block.  A  group  of  five  short  wire  segments,  cut 
from  the  same  size  of  wire,  are  placed  near  the  tip  of  each  wire.  These  groups  of  short 
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wire  segments  are  intended  for  simulating  microcalcifications.  They  were  also  placed  in 
overlap  with  the  anthropomorphic  phantom  to  simulate  microcalcifications  in  the  breast. 
This  phantom  allows  subtle  difference  between  different  imaging  systems  or  techniques 
to  be  detected. 

Image  Acquisition 

The  three  phantoms  were  first  imaged  with  a  conventional  screen-film 
combination  (Kodak  Min-R  screen-film  combination)  in  the  photo-timed  mode.  After 
ensuring  that  the  resulting  film  has  adequate  density,  a  cassette  containing  a  CR  screen 
was  exposed  with  exactly  the  same  technique.  All  images  were  taken  at  25  kVp  with  the 
reciprocating  grid  in  place.  The  mAs  used  were  144,  22  and  104  for  the  tissue  equivalent 
breast  phantom,  perception  phantom  and  anthropomorphic  breast  phantom  respectively. 
Two  different  types  of  CR  screens,  HR  and  ST,  were  used  for  exposure.  The  exposed 
screens  were  scanned  in  the  high  resolution  mode  using  a  pixel  size  of  43  pm  and  a  laser 
beam  spot  size  of  30  or  50pm  (full  width  at  half-maximum). 

Phantom  Images 

Films  5  shows  a  screen-film  image  of  the  perception  phantom  acquired  with  a 
mammographic  unit  at  25  kVp  and  22  mAs.  Films  6, 7  and  8  show  CR  images  of  the 
same  phantom  acquired  with  three  different  techniques:  (Film  6)  exposed  with  a  ST 
screen  and  read  out  with  30pm  laser  beam  size,  (Film  7)  exposed  with  a  HR  screen  and 
read  out  with  30pm  laser  beam  size  and  (Film  8)  exposed  with  a  HR  screen  and  read  out 
with  50pm  laser  beam  size.  All  images  were  acquired  at  25  kVp  and  22  mAs  and  read 
out  with  a  pixel  size  of  43  pm. 

The  wires  segments  have  a  sharper  appearance  in  the  SF  image.  In  the  screen- 
film  image  (Film  5),  the  5  mil  or  thicker  (125pm)  wire  segments  are  well  resolved  while 
the  4  mil  (100pm)  wire  segments  are  only  marginally  resolved.  In  the  optimized  CR 
image  (Film  7),  all  wires  with  a  thickness  of  5mil  or  thicker  are  resolved.  The  long  wires 
are  easier  to  see.  All  but  the  2  mil  (50  pm)  long  wires  (oriented  in  radial  direction)  are 
resolved.  The  screen-film  image  shows  the  wire  segments  with  a  sharper  appearance 
which  is  partially  due  to  its  smaller  (by  a  factor  of  2)  format.  This  helps  resolve  the  4  mil 
(100pm)  wire  segments.  However,  it  does  not  resolve  finer  segments  (2,  3  mil)  or  wires 
(2  mil).  This  shows  that  at  typical  clinical  exposure  level  (as  determined  by  the 
automatic  exposure  control),  the  improved  storage  phosphor  imaging  system  approaches 
the  conventional  screen-film  combinations  in  quality  as  far  as  detection  of 
microcalcifications  is  concerned.  The  use  of  larger  beam  spot  size  or  thicker  screen  does 
not  seem  to  affect  the  detection  of  the  5mil  or  thicker  wires  in  the  CR  images. 

Films  9  shows  a  screen-film  image  of  the  tissue  equivalent  breast  phantom 
acquired  with  a  mammographic  unit  at  25  kVp  and  144  mAs.  Films  10,  1 1  and  12  show 
CR  images  of  the  same  phantom  acquired  with  three  different  techniques:  (Film  2) 
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exposed  with  a  ST  screen  and  read  out  with  30pm  laser  beam  size,  (Film  3)  exposed  with 
a  HR  screen  and  read  out  with  30pm  laser  beam  size  and  (Film  4)  exposed  with  a  HR 
screen  and  read  out  with  50pm  laser  beam  size.  All  images  were  acquired  at  25  kVp  and 
144  mAs  and  read  out  with  a  pixel  size  of  43pm. 

Probably  because  of  their  spherical  shape,  the  simulated  microcalcifications  in  the 
tissue  equivalent  breast  phantom  are  more  difficult  to  detect  than  the  aluminum  wire 
segments  in  the  perception  phantom.  Both  SF  and  CR  images  can  resolve  200pm 
microcalcifications.  Neither  can  resolve  the  160pm  microcalcifications.  Because  the  CR 
image  can  be  enhanced  for  better  display  contrast,  the  SP  image  actually  shows  the 
circular  simulated  masses  better  than  the  screen-film  image.  However,  it  doesn’t  resolve 
smaller  masses  than  the  screen-film  image.  The  screen-film  image  shows  the 
calcifications,  masses  and  fabrils  with  a  sharper  appearance  which  is  partially  due  to  its 
smaller  (by  a  factor  of  2)  format.  However,  it  does  not  resolve  more  microcalcifications 
than  what  can  be  resolved  in  the  CR  image.  It  seems  that  the  size  differences  of  the 
simulated  microcalcifications  and  masses  in  this  phantom  are  too  coarse  to  differentiate 
the  performance  differences  between  the  SF  and  various  CR  techniques. 

Films  13  shows  a  screen-film  image  of  an  anthropomorphic  breast  phantom 
acquired  with  a  mammographic  unit  at  25  kVp  and  104  mAs.  Films  2,  3  and  4  show  CR 
images  of  the  same  phantom  acquired  with  three  different  techniques:  (Film  14)  exposed 
with  a  ST  screen  and  read  out  with  30pm  laser  beam  size,  (Film  15)  exposed  with  a  HR 
screen  and  read  out  with  30pm  laser  beam  size  and  (Film  16)  exposed  with  a  HR  screen 
and  read  out  with  50pm  laser  beam  size.  All  images  were  acquired  at  25  kVp  and  104 
mAs  and  read  out  with  a  pixel  size  of  43  pm. 

With  the  SF  imageg,  5  mil  or  thicker  wires  can  be  detected  while  with  the  best  CR 
image  (Film  15)  only  6  mil  and  thicker  wires  can  be  detected.  Notice  that  with  the  CIRS 
breast  phantom  or  the  perception  phantom,  the  contrast  objects  are  detected  over  uniform 
background  and  the  SP  and  SF  are  comparable  to  each  other  in  performance.  The  sharper 
appearance  of  the  SF  image  has  an  advantage  here  probably  because  the  sharper 
appearance  of  the  SF  image  help  differentiate  the  simulated  microcalcifications  from  the 
overlying  background  structures.  The  use  of  thicker  screen  or  larger  beam  spot  size  result 
in  more  blurried  look  of  the  images,  making  it  slightly  more  difficult  to  see  the  simulated 
microcalcifications. 

There  are  several  lessons  learned  from  this  comparison  study  based  on  visual 
comparison  of  the  phantom  images.  First  of  all,  commercial  phantoms  are  inadequate  for 
comparing  the  CR  with  the  SF  techniques.  The  problems  seem  to  reside  on  the  fact  that 
the  size  differences  in  the  simulated  microcalcifications  and  masses  are  too  coarse  to  be 
used  to  depict  any  subtle  difference  in  performance.  The  use  of  aluminum  wire  segments 
with  a  range  of  closely  spaced  sizes  help  depict  the  difference.  However,  a  far  better 
comparison  technique  should  be  the  use  of  ROC  studies. 
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COMPARISON  STUDY-  ROC  STUDIES 


The  ROC  studies  provide  a  means  to  quantitatively  compare  two  imaging 
techniques.  Two  ROC  studies  were  conducted  to  compare  the  improved  CR  imaging 
technique  with  the  conventional  screen-film  technique.  The  studies  were  based  on 
phantom  images  with  simulated  microcalcifications  and  masses.  To  present  more 
realistic  imaging  conditions,  these  objects  were  superimposed  with  structures  mimicking 
tissue  structures  of  breasts.  Methods  and  results  for  these  are  described  and  discussed  in 
the  following  sections. 

Phantom  Design 

The  main  body  of  the  phantom  is  supposed  to  provide  a  background  structure  over 
which  microcalcifications  and  masses  will  be  simulated  and  superimposed.  We  have 
investigated  several  approaches  to  construct  the  main  body  of  the  phantom.  The  final 
design  is  a  large  1"  thick  piece  of  beef  housed  in  a  specifically  designed  plastic  tank. 

This  design  is  intended  to  simulate  a  well  compressed  breast.  The  meat  produces  a 
structural  pattern  different  form  real  mammograms  but  still  provides  adequate  challenge 
to  the  tasks  of  detecting  simulated  microcalcifications  or  masses.  The  meat  can  be  easily 
replaced  to  produce  a  different  structural  background  so  that  the  reader  has  no 
opportunity  to  become  "trained"  for  a  certain  pattern.  The  tank  is  filled  with  pure  alcohol 
to  preserve  the  meat  for  longer  time.  The  thiekness  of  the  tank  is  minimized  to  about 
0.5".  Thus,  the  phantom  has  a  total  thickness  of  1 .5". 

Mircocalcifications  (MC)  are  simulated  by  a  group  of  5  or  6  aluminum  wire 
segments  (spread  over  an  area  of  2-3  mm  wide)  placed  on  a  a  1 .5"xl.5"  piece  of  plastic 
film  (cut  from  acetone  transparency).  The  MCs  are  simulated  with  two  different  sizes:  5 
or  6  mils(125  or  150pm).  Other  sizes  weren't  used  because  visual  tests  have  indicated 
that  MCs  larger  than  6  mil  are  too  easy  to  see  while  those  smaller  than  5  mil  can  never  be 
seen  when  overlapping  with  some  tissue  structures.  For  each  size,  many  film  inserts  were 
made  with  the  MCs  located  at  various  locations.  These  films  can  also  be  rotated  by 
multiples  of  90  degrees  to  further  randomize  the  location  of  the  MC  groups.  Films 
without  any  MC  were  also  made  to  simulate  normal  cases.  A  film  holder  is  constructed 
of  a  3/8"  thick  piece  of  Lucite  with  a  3x4  matrix  of  1.5"xl.5"  holes  cut  out.  During 
image  acquisition,  the  holder,  with  12  selected  film  inserts,  is  placed  underneath  the  meat 
phantom  for  exposure.  Figure  13  shows  the  layout  of  these  cutouts.  Films  17  and  18 
show  a  CR  and  a  screen-film  image  of  the  meat  phantom  with  simulated  MCs. 

Statistical  analysis  in  the  proposal  has  indicated  that  240  observations  (readings) 
are  required  to  yield  reasonable  significance  of  the  ROC  study.  Our  phantom  design 
provides  12  "windows".  Within  each  of  these  windows,  an  observation  or  reading  ean  be 
performed  to  determine  whether  and  where  there  are  MCs.  Thus,  each  image  will 
provide  12  independent  readings.  In  our  studies,  a  total  of  24  images  were  acquired, 
providing  a  total  of  288  different  detection  tasks  for  each  reader.  144  film  inserts  were 


31 


None 


1 


2 


3 


4  5  6 


Figure  OLayout  of  the  film  holder  for  the  meat  phantom  used  to  acquire  images  for  the 
ROC  study.  A  matrix  of  3x4  1.5"xl.5"  film  inserts  can  be  placed  on  the  holder 
to  provide  12  independent  reading  tasks.  Examples  of  the  microcalcification 
placement  are  shown  in  the  windows  #2,3,5,7,8,9,10,1 1,12. 
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generated  without  any  microcalcifications  and  used  as  control.  72  film  inserts  were 
produced  with  6  mil  MC’s.  Another  72  film  inserts  were  produced  with  5  mil  MC’s. 
Using  a  computer  program,  these  inserts  were  randomly  selected  and  placed  during  the 
acquisition  of  each  image. 

The  masses  are  simulated  by  a  circular  piece  of  film..  The  edges  of  the  film 
pieces  were  ground  to  obtain  smooth  border.  The  variation  of  the  grinding  procedure 
produces  mass-like  objects  with  variable  and  marginal  visibility.  A  total  of  24  images 
were  acquired,  each  containing  3x4  1.5”xl.5”  small  regions  for  12  independent  detection 
tasks.  Among  a  total  of  288  image  regions,  half  of  them  have  no  simulated  mass  and  half 
of  them  have  one  mass  placed  at  the  center  of  the  image  or  at  center  of  one  of  the  four 
quadrants.  The  distribution  and  placement  of  the  massess  in  the  images  were  randomized 
with  the  help  of  a  computer  program.  Films  19  and  20  show  a  CR  and  a  screen-film 
image  of  the  meat  phantom  with  simulated  masses. 


Image  Acquisition 

A  clinical  mammographic  unit  (LoRAD  M-III)  is  used  for  image  acquisition.  The 
x-ray  tube  uses  a  molybdenum  target.  A  5:1  reciprocating  antiscatter  grid  is  used.  With 
the  phantom  in  place,  a  screen-film  image  is  first  taken  using  Kodak  Min-R  screen  and 
Min-R  E  film  in  the  autotimed  mode.  The  kVp  is  kept  fixed  at  25  kVp  and  the  mAs  is 
controlled  by  photo-timing.  Following  the  exposure,  the  mAs  will  be  recorded.  The 
exposed  film  will  then  be  developed  and  checked  for  proper  exposure.  If  the  film 
exposure  is  adequate,  the  storage  phosphor  cassette  will  be  inserted  and  exposed  using 
identical  x-ray  techniques.  The  exposed  storage  phosphor  screen  will  then  be  scanned 
with  a  pixel  size  of  43pm  to  generate  a  4096x5000x12  bit  digital  image.  The  image  will 
first  be  printed  with  a  80pm  pixel  size  using  a  Kodak  Ektascan  laser  printer.  The  film 
will  be  examined  and  checked  for  any  artifacts  or  phantom  position.  If  acceptable,  the 
phantom  will  be  fitted  with  a  different  set  of  film  inserts  and  imaged  again  following  the 
same  procedure. 

Data  Collection 


The  CR  images  are  printed  with  a  pixel  size  of  50  micron  to  produce  films  of  a 
size  (8”xl0”)  similar  to  that  of  mammograms  (18cm  x  24cm).  Because  we  do  not  have  a 
printer  capable  of  printing  with  50pm  pixel  size,  the  CR  image  data  were  sent  to  Eastman 
Kodak  Company  for  printing.  Within  each  ROC  study,  either  microcalcification  or  mass 
detection  study,  a  total  of  24  x  12  =  288  sub-images  are  available  for  reading  and  scoring 
for  each  imaging  technique.  Five  board  certified  mammographers  were  asked  to  read 
and  score  the  288  sub-images  acquired  with  the  CR  technique  and  another  288  sub¬ 
images  acquired  with  the  SF  technique..  A  computer  program  was  written  to  manage  the 
reading  study  and  to  record  the  scores.  The  scores,  ranging  from  0  to  1 00,  represent  the 
confidency  level  for  the  microcalcifications  or  mass  to  be  present.  If  a  score  of  50  or 
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higher  was  entered,  the  reader  was  also  asked  to  enter  the  perceived  location  of  the 
microcalcification  group  or  mass.  This  location  information  was  recorded  along  with  the 
confidence  levels  for  future  analysis. 

ROC  Data  Analysis 

The  results  fi'om  the  reading  sessions  were  analyzed  using  the  ROC  software 
developed  by  Dr.  Metz  at  the  University  of  Chicago.  Notice  that  in  ROC  analysis  the 
location  information  is  ignored.  The  location  information  will  be  used  in  the  LROC 
analysis  which  would  lower  the  scores  if  the  locations  entered  are  incorrect  ones.  It 
would  be  interesting  to  compare  the  LROC  analysis  with  the  ROC  analysis  to  see  if 
incorrect  locations  has  any  impact  on  the  final  outcome.  In  a  small  scale  preliminary 
study,  we  have  found  that  the  locations  entered  by  the  readers  are  mostly  correct  ones  and 
therefore  there  is  no  need  to  use  the  location  information  to  perform  the  LROC  analysis. 

ROC  curves  were  formed  by  plotting  the  true  positive  fraction  against  the  false 
positive  fraction  for  each  confidence  level.  A  ROC  curve  was  computed  for  each 
imaging  technique  and  for  each  reader.  The  area  underneath  the  ROC  curve  was  also 
computed  for  quantitatively  comparing  the  performances  of  the  CR  and  SF  techniques. 
The  results  indicate  that  for  detection  of  the  microcalcifications,  the  conventional  screen- 
film  mammography  technique  outperforms  our  high  resolution  CR  system.  However,  for 
detection  of  the  masses,  our  high  resolution  CR  system  outperforms  the  screen-film 
system. 


Figures  14  and  15  show  two  examples  of  the  ROC  curves  obtained  for  detection 
of  the  microcalcifications.  In  each  figure,  the  ROC  curve  for  CR  is  plotted  along  with 
that  for  the  SF  technique  for  comparison.  Each  figure  corresponds  to  a  different  reader. 
The  areas  under  the  ROC  curves  are  0.9060  for  SF  versus  0.7259  for  CR  for  Reader  1 
(Figure  14)  and  0.8368  for  SF  versus  0.7366  for  CR  for  Reader  2  (Figure  15).  The  results 
clearly  shows  that  despite  variation  among  readers  the  SF  images  perform  better  than  the 
CR  images  in  detecting  microcalcifications. 

A  similar  ROC  study  for  detection  of  simulated  masses  has  also  been  conducted. 
Data  from  a  preliminary  ROC  study  were  similarly  analyzed.  Figures  16  shows  an 
example  of  the  ROC  curves  for  the  detection  of  masses.  The  areas  under  the  ROC  curves 
are  0.9655  for  the  CR  technique  versus  0.8270  for  SF  technique.  The  results  show  that 
for  detction  of  masses  the  CR  images  perform  significantly  better  than  the  SF  images. 
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Figure  14  ROC  curves  for  dectecting  the  simulated  microcalcifications:  computed 
radiography  versus  screen-film  mammography  techniques  (Example  1). 
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Figure  15  ROC  curves  for  dectecting  the  simulated  microcalcifications:  computed 
radiography  versus  screen-film  mammography  techniques  (Example  2). 
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Figure  16  ROC  curves  for  dectecting  the  simulated  microcalcifications:  computed 
radiography  versus  screen-film  mammography  techniques  (Example  1). 
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COMPRESSION  OF  IMAGE  DATA 

Two  methods  were  used  to  develop  and  implement  algorithms  for  compression  of 
digital  mammographic  image  data.  The  first  method  is  a  lossless  scheme  which  employs 
image  segmentation  technique  to  eliminate  data  outside  the  breast  area.  The  second 
method  uses  JPEG  image  data  compression  algorithm  to  reduce  image  data  in  a  lossy 
fashion.  Although  this  algorithm  has  been  investigated  for  image  data  compression  in 

other  imaging  applications^^,  it  must  be  customized  for  application  to  mammographic 
images.  This  algorithm  also  results  in  irreversible  change  of  image  data.  Therefore, 
comparison  studies  were  performed  to  determine  the  parameters  used  in  the  algorithm. 
Upon  development  and  optimization  of  the  algorithm,  more  vigorous  comparison  studies 
will  be  performed  to  validate  the  algorithm  by  ensuring  that  no  significant  loss  of  image 
quality  has  been  incurred. 

The  data  compression  procedure  developed  so  far  can  be  divided  into  two  parts: 
preprocessing  and  Joint  Photographic  Experts  Group  (JPEG)  compression.  They  are 
summarized  in  this  section.  A  detailed  description  and  discussion  of  the  methods  are 
presented  in  the  manuscript  reprint  attached  as  Appendix  V. 

Preprocessing 

Before  applying  the  JPEG  compression  algorithm,  the  image  data  is  preprocessed 
to  improve  its  compressibility.  This  preprocessing  operation  is  divided  into  two  phases; 
first  the  image  is  segmented  and  cropped  and  the  background  pixel  values  are  modified  to 
minimize  the  storage  requirements  for  the  background.  We  begin  the  segmentation 
process  by  automatically  determining  a  threshold  value  from  the  histogram  of  the 
digitized  source  image  that  indicates  the  transition  between  background  values  and  tissue 
values.  This  threshold  value  is  used  to  classify  pixels  in  the  image  as  either  tissue  or 
background  so  as  to  produce  a  two-valued  mask  corresponding  to  the  tissue-background 
separation. 

The  second  phase  of  preprocessing  involves  the  application  of  a  noise-removal 
filter  to  the  tissue  area  of  the  image.  The  particular  filter  tested  here  is  a  nonlinear  order 
filter  that  was  designed  to  modify  mainly  pixels  with  large  noise  components  while 
leaving  the  majority  of  pixels  unchanged. 

JPEG  Compression 

The  JPEG  compression  algorithm  used  in  this  study  is  a  software  package  written 
and  tested  in  our  laboratory  that  implements  the  12-bit  version  of  the  extended  JPEG 
standard.  For  this  study,  Huffman  encoding  was  used  to  encode  the  quantized 
coefficients  and  all  compression  ratios  were  based  on  this.  The  Huffman  tables  were 
derived  from  the  statistics  of  each  individual  quantized  image  rather  than  on  the  ensemble 


38 


C.  Shaw,  DAMD  17-93-J-3009 


statistics.  Signaling  information  was  not  included  in  the  code  size  when  calculating 
compression  ratios. 

Evaluation 


The  compression  technique  is  evaluated  by  its  speed,  compression  ratio  and 
ability  to  preserve  image  information.  Measurement  of  the  compression  and 
decompression  speed  is  straightforward.  Ther  results,  however,  are  dependent  upon 
individual  images.  Therefore,  we  evaluated  and  compared  the  average  speeds  for  the  8 
images  used  in  our  study. 

The  purpose  of  preprocessing  is  to  increase  the  information  content  in  the  image. 
Thus,  we  measured  the  average  number  of  bits  for  each  8x8  pixel  block  of  tissue  area 
before  and  after  segmentation.  An  improvement  factor  is  computed  as  the  ratio  of 
average  number  of  bits  before  the  segmentation  to  that  after  the  segmentation.  We  also 
measured  the  percentage  of  image  area  with  breast  tissue  signals. 

The  compression  ratio  and  ability  for  information  preservation  are  not 
independent  of  each  other.  It  is  reasonable  to  expect  better  information  preservation 
when  the  compression  ratio  is  low  and  vice  versa.  We  used  a  just  noticeable  difference 
(JND)  study  to  quantify  the  information  preservation  as  a  function  of  the  compression 
ratio.  From  eight  original  mammograms,  which  had  each  been  compressed  at  five 
compression  ratios,  a  set  of  40  randomized  pairs  of  images,  each  consisting  of  one  laser 
printed  non-compressed  image  and  one  compressed  image  were  assembled.  20  readers 
were  asked  to  view  the  pairs  side  by  side  and  forced  to  choose  the  compressed  image. 
The  readers  include  experienced  radiologists,  physicists  and  engineers  specializing  in 
medical  imaging.  Following  the  study,  the  results  are  sorted  according  the  compression 
ratio  and  the  percentage  of  corrected  selected  pairs  is  plotted  as  a  function  of  the 
compression  ratio  for  a  group  of  "expert  readers"  and  a  group  of  "average  readers". 

We  will  summarize  our  results  from  the  data  compression  study  in  regard  to  three 
topics:  speed  of  our  compression  system,  performance  of  preprocessing  and  performance 
of  the  JPEG  compression.  A  manuscript  is  attached  as  Appendix  V  for  a  more  detailed 
description  and  discussion  of  the  results. 

Speed 


The  speed  of  our  compression  system  is  somewhat  image  dependent.  It  averages 
about  20  s/Mbyte  to  preprocess  and  compress  source  data.  Decompression  is  faster  at  1 5 
s/Mbyte  because  there  is  no  preprocessing  involved.  These  speeds  may  not  be  acceptable 
for  clinical  use  with  a  picture  archival  and  communication  system  (PACS).  However, 
they  can  be  easily  improved  by  using  dedicated  and  optimized  hardware  as  much  as 
possible. 
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Preprocessing 

A  summary  of  the  effect  of  image  segmentation  on  the  number  of  bits  needed  to 
represent  tissue  blocks  is  shown  in  Table  1  in  which  the  percentage  of  image  area 
containing  tissue  signals,  bits  per  tissue  block  before  and  after  segmentation,  and  the 
improvement  factor  are  listed  for  each  of  the  8  original  mammograms.  From  these 
results,  the  information  content  of  the  image  has  been  improved  by  a  factor  of  2.7  on  the 
average  despite  the  fact  that  the  background  only  represented  36%  of  the  image  area. 

JPEG  Compression 

In  Figure  17,  the  performance  of  image  processing  specialists  and  nonspecialists, 
measured  by  the  ratio  of  correct  identification  of  compressed  images,  is  plotted  as  a 
function  of  the  compression  ratio.  These  results  verify  our  expectation  that  the  higher  the 
compression  ratio  is,  the  more  information  change  there  is  and  therefore  the  easier  one 
can  identify  the  compressed  images.  However,  the  impact  of  data  compression  on 
diagnostic  accuracy  must  be  evaluated  by  a  ROC  study. 
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Table  1:  The  Effect  of  Image  Segmentation  on  Tissue  Block 
Representation  {bits  per  block) 


Image 

%  of  Area 
Representing 
Tissue 

Bits  per 
Tissue  Block 
(no 

preprocessing) 

Bits  per 
Tissue  Block 
(after 

segmentation) 

Improvement 

Factor 

1 

51.1 

8.6 

30.1 

3.5 

2 

52.9 

8.5 

27.0 

3.2 

3 

70.5 

9.1 

32.3 

3.5 

4 

71.7 

10.0 

28.0 

2.8 

5 

72.4 

9.4 

23.6 

2.5 

6 

56.8 

7.6 

26.5 

3.5 

7 

70.7 

24.2 

27.4 

1.1 

8 

67.0 

18.0 

22.2 

1.2 

Fig  17  Performance  of  image-processing  specialists  and 
nonspecialists  for  the  task  of  identifying  compressed  versus 
noncompressed  images.  41 
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CONCLUSIONS 


To  summary  our  accomplishment  in  the  grant  period: 

1 .  We  have  developed  and  implemented  a  technique  to  measure  and  monitor  the  laser 
beam  spot  size  (Half  Maximum  Width).  Furthermore,  the  operation  was  automated 
and  sped  up  so  that  the  laser  beam  profile  and  spot  size  can  be  monitored  during 
focusing  or  adjustment  of  the  beam  size.  With  the  help  of  this  technique,  we  have 
successfully  modified  the  optics  of  our  storage  phosphor  image  reader  and  reduced  the 
laser  beam  spot  size  to  approximately  30|am  (Half  Maximum  Width). 

2.  We  have  measured  and  studied  the  characteristics  of  of  our  CR  system.  A  Fourier 
transform  technique  was  developed  to  measure  the  MTF  of  digital  imaging  system 
with  images  of  commercial  resolution  bar  patterns.  Using  this  technique,  the  MTF’s 
of  our  system  were  measured  and  compared  with  those  of  a  typical  commercial 
system.  Effects  of  reduced  laser  beam  spot  and  CR  screens  were  also  studied.  A 
quality  control  protocal  has  been  developed  and  used  to  check  the  performance  of  our 
CR  system  during  the  course  of  our  project. 

3.  We  have  developed  a  perception  phantom  for  evaluating  and  monitoring  the  low 
contrast  performance  of  mammographic  imaging  systems.  Using  this  phantom  and 
two  other  phantoms,  a  series  of  images  were  acquired  and  used  to  study  the  effects  of 
reduced  laser  beam  spot  and  screen  selection  in  CR  mammography  and  to  compare  the 
CR  with  SF  systems  for  mammographic  imaging. 

4.  We  have  designed  and  constructed  phantoms  for  the  ROC  studies.  Using  these 
phantoms,  we  have  acquired  both  CR  and  SF  images  to  simulate  the  task  for  detecting 
microcalcifications  and  the  task  for  detecting  masses.  The  CR  images  were  printed 
with  a  pixel  size  of  SOpm  for  the  reading  studies.  Software  for  acquiring  and 
analyzing  the  reading  studies  were  developed.  Reading  experiments  have  been 
conducted.  ROC  curves  have  been  obtained  to  compare  the  CR  with  the  SF 
techniques. 

5.  We  have  developed  a  JPEG  based  algorithm  for  compressing  meimmographic  image 
data.  A  Just  Noticeable  Difference  (JND)  study  was  conducted  to  quantify  the 
information  preservation  as  a  function  of  the  compression  ratio. 

The  above  studies  have  resulted  the  following  conclusions: 

1 .  We  have  demonstrated  it  is  feasible  to  configure  a  CR  system  for  high  resolution 
imaging  applications  like  mammography.  It  is  possible  to  modify  the  laser  scanning 
optics  to  reduce  the  laser  beam  spot  size  down  to  40p,m  or  even  smaller. 
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2.  The  reduction  of  the  laser  beam  spot  size  and  use  of  a  high  resolution  CR  screen 
result  in  a  significant  improvement  of  the  spatial  resolution  capability  of  the  storage 
phosphor  imaging  system.  This  improvement  allows  us  to  visually  resolve  10  Ips/mm 
bar  patterns  at  25  kVp  on  a  mammographic  unit.  The  MTF’s  of  the  improved  CR 
system  significantly  exceed  those  of  a  typical  commercial  system.  The  greater 
detection  efficiency  from  using  a  thicker  screen  (ST  screen)  is  more  than  offset  by  the 
lossofMTF. 

3.  Commercial  phantoms  alone  can  not  differentiate  the  subtle  performance  differences 
between  the  CR  and  SF  techniques.  The  use  of  aluminum  wires  with  finely  spaced 
sizes  help  differentiate  the  performance  difference  between  the  CR  and  SF 
techniques.  Comparison  with  phantom  images  show  that  the  performance  of  our  CR 
system  approach  but  is  slightly  lower  than  that  of  a  conventional  SF  mammography 
system.  However,  from  the  perspective  of  visual  appearance,  the  difference  is  rather 
subtle. 

4.  Our  ROC  studies  indicate  that  for  the  detection  of  microcalcifications,  the  SF  images 
perform  better  than  the  CR  images.  For  detction  of  masses,  the  CR  images  seem  to 
perform  better  than  the  SF  images. 

5.  The  results  of  our  data  compression  studies  indicate  that  the  detectable  information 
change  increases  with  compression  ratio.  However,  the  visual  detection  of 
information  change  may  not  be  directly  related  to  diagnosis.  Direct  impact  of  data 
compression  on  diagnostic  accuracy  must  be  evaluated  by  a  ROC  study  which  is 
beyond  the  scope  of  our  image  compression  study. 

The  implications  of  these  conclusion  are  as  follows: 

While  our  CR  system  has  been  significantly  improved  by  reducing  the  laser  beam 
spot  size,  it  could  benefit  from  further  improvement  of  spatial  resolution  for  the  task  of 
detecting  small  microcalcifications  (less  than  150pm  in  dimension).  However,  detection 
of  such  small  microcalcifications  is  rather  infrequent  in  conventional  mammography. 
Therefore,  its  role  in  screening  brast  cancer  is  questionable.  If  detection  of 
microcalcifications  with  a  dimension  of  smaller  than  150pm  is  essential  to  diagnosis, 
then  further  improvement  of  the  spatial  resolution  quality  of  CR  imaging  is  called  for.  As 
the  laser  beam  spot  size  has  been  greatly  reduced  to  an  extent  that  the  MTF  degradation  is 
rather  insignificant,  further  increase  of  the  MTF  should  be  achieved  by  improving  the 
design  of  CR  image  plates  used  for  exposure  rather  than  further  reduction  of  the  laser 
beam  spot  size.  This  requies  the  development  of  either  a  more  efficient  screen  material  to 
reduce  the  screen  thickness  or  a  structured  CR  screen  to  eliminate  or  reduce  lateral 
difussion  of  the  laser  beam. 
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SPECIFICATIONS  FOR  THE  RESOLUTION  TEST  PATTERN 
(Nuclear  Associates,  Model  07-553) 


AI-1 


AI-2 


APPENDIX  II 


SPECIFICATIONS  FOR  THE  CIRS  BREAST  PHANTOM 
(Computerized  Imaging  Reference  Systems,  Inc.) 


AII-l 


Map  of 

Embedded  Details  in  Phantoms, 


•  Stepwedge 

6  mm  glandular  equivalent/5  mm  fat 
equivalent 


•  Masses 

6  simulated  tumor  masses  ranging  from 
6  mm  diameter  to  1  mm  diameter 


•  Microcalcifications  (CaCOJ  •  Fibrils 

1 2  groupings  with  largest  at  .35  mm  3  fibre  spiculations 

diameter  and  smallest  at  .12  mm 

diameter.  This  range  of  sizes  has  been 

shown  in  practice  to  be  sensitive  to 

system  calibration  and  performance. 


rJTCTrB  COMPUTERIZED  IMAGING 
REFERENCE  SYSTEMS,  INC. 

2488  Almedo  Avenue.  Norfolk,  Virginia  23513 
Tel:  (804)855-2765 
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APPENDIX  III 


SPECIFICATIONS  FOR  THE  ANTHROPOMORPHIC  BREAST  PHANTOM 
(Model  165,  Gammex-Radiation  Measurement  Inc.) 


AIII-1 


Specifications 


Breast  Tissue  Equivalent 
Plastic  Casting 


High  Resolution  Layer 


Phantom  Case 


Overall  Dimensions 


Tissue  Equivalent  Plastic, 

RMI  Model  454 

50%  adipose,  50%  glandular, 

(BR-12  equivalent) 

Radiographic  film  containing 
a  stable  mercury-silver 
amalgam. 

9  step  gray  level  stepwedge. 

Resolution  Bar  Pattern: 

5-25  Ip/mm 

Case  Material:  Acrylic 

Top  Thickness:  0.215  in 
Bottom  Thickness:  0.175  in 

19.6  X  11.7  X  6.1  cm 
(7J  X  4.7  X  2.4  in) 

0.7  kg  (1.61b) 


Note:  Due  to  our  philosophy  of  continuous  product 
improvement,  these  specifications  may 
change  without  notice. 
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Film  1 


Film  2 


Film  3 


Film  4 


Film  5 


Film  6 


Film  7 


Film  8 


Film  9 


Film  10 


APPENDIX  IV 
Captions  for  Attached  Films 


Image  of  a  resolution  bar  pattern  obtained  with  a  screen-film  mammography 
system.  Exposures  were  made  at  25kVp,  6mAs  and  a  source-to-image  distance 
of  65cm. 

CR  image  of  a  resolution  bar  pattern  obtained  with  a  ST  screen,  a  laser  beam 
diameter  of  30pm  and  a  pixel  size  of  43pm  were  used.  Exposures  were  made  at 
25kVp,  lOmAs  and  a  source-to-image  distance  of  65cm. 

CR  image  of  a  resolution  bar  pattern  obtained  with  a  HR  screen,  a  laser  beam 
diameter  of  30pm  and  a  pixel  size  of  43pm  were  used.  Exposures  were  made  at 
25kVp,  lOmAs  and  a  source-to-image  distance  of  65cm. 

CR  image  of  a  resolution  bar  pattern  obtained  with  a  HR  screen,  a  laser  beam 
diameter  of  50pm  and  a  pixel  size  of  43pm  were  used.  Exposures  were  made  at 
25kVp,  lOmAs  and  a  source-to-image  distance  of  65cm. 

Image  of  a  perception  phantom  obtained  with  a  screen-film  mammography 
system.  Exposures  were  made  at  25kVp,  22mAs  and  a  source-to-image 
distance  of  65  cm. 

CR  image  of  a  perception  phantom  obtained  with  a  ST  screen,  a  laser  beam 
diameter  of  30pm  and  a  pixel  size  of  43pm  were  used.  Exposures  were  made  at 
25kVp,  22mAs  and  a  source-to-image  distance  of  65cm. 

CR  image  of  a  perception  phantom  obtained  with  a  HR  screen,  a  laser  beam 
diameter  of  30pm  and  a  pixel  size  of  43pm  were  used.  Exposures  were  made  at 
25kVp,  22niAs  and  a  source-to-image  distance  of  65cm. 

CR  image  of  a  perception  phantom  obtained  with  a  HR  screen,  a  laser  beam 
diameter  of  50pm  and  a  pixel  size  of  43  pm  were  used.  Exposures  were  made  at 
25kVp,  22mAs  and  a  source-to-image  distance  of  65cm. 

Image  of  a  tissue  equivalent  breast  phantom  obtained  with  a  screen-film 
mammography  system.  Exposures  were  made  at  25kVp,  144mAs  and  a  source- 
to-image  distance  of  65  cm. 

CR  image  of  a  tissue  equivalent  breast  phantom  obtained  with  a  ST  screen,  a 
laser  beam  diameter  of  30pm  and  a  pixel  size  of  43  pm  were  used.  Exposures 
were  made  at  25kVp,  144mAs  and  a  source-to-image  distance  of  65cm. 
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Film  1 1  CR  image  of  a  tissue  equivalent  breast  phantom  obtained  with  a  HR  screen,  a 
laser  beam  diameter  of  30pm  and  a  pixel  size  of  43pm  were  used.  Exposures 
were  made  at  25kVp,  144mAs  and  a  source-to-image  distance  of  65cm. 

Film  12  CR  image  of  a  tissue  equivalent  breast  phantom  obtained  with  a  HR  screen,  a 
laser  beam  diameter  of  50pm  and  a  pixel  size  of  43  pm  were  used.  Exposures 
were  made  at  25kVp,  144mAs  and  a  source-to-image  distance  of  65cm. 

Film  1 3  Image  of  an  anthropomorphic  breast  phantom  obtained  with  a  screen-film 

mammography  system.  Exposures  were  made  at  25kVp,  104mAs  and  a  source- 
to-image  distance  of  65cm. 

Film  14  CR  image  of  an  anthropomorphic  breast  phantom  obtained  with  a  ST  screen,  a 
laser  beam  diameter  of  30pm  and  a  pixel  size  of  43pm  were  used.  Exposures 
were  made  at  25kVp,  104mAs  and  a  source-to-image  distance  of  65cm. 

Film  1 5  CR  image  of  an  anthropomorphic  breast  phantom  obtained  with  a  HR  screen,  a 
laser  beam  diameter  of  30pm  and  a  pixel  size  of  43  pm  were  used.  Exposures 
were  made  at  25kVp,  104mAs  and  a  source-to-image  distance  of  65cm. 

Film  16  CR  image  of  an  anthropomorphic  breast  phantom  obtained  with  a  HR  screen,  a 
laser  beam  diameter  of  50pm  and  a  pixel  size  of  43pm  were  used.  Exposures 
were  made  at  25kVp,  104mAs  and  a  source-to-image  distance  of  65cm. 

Film  17  CR  image  of  a  meat  phantom  with  simulated  microcalcifications.  Exposures 
were  made  at  25kVp,  22mAs  and  a  source-to-image  distance  of  65cm.  5mil 
microcalcifications  are  present  in  box  #1, 2,  5,  7  and  6mil  microcalcifications 
were  simulated  in  box  #8, 11.  The  boxes  are  numbered  as  shown  in  Figure  13. 

Film  18  SF  image  of  a  meat  phantom  with  simulated  microcalcifications.  Exposures 
were  made  at  25kVp,  22mAs  and  a  source-to-image  distance  of  65cm.  5mil 
microcalcifications  are  present  in  box  #1,  2,  5,  7  and  6mil  microcalcifications 
were  simulated  in  box  #8, 11.  The  boxes  are  numbered  as  shown  in  Figure  1 3 . 

Film  19  CR  image  of  a  meat  phantom  with  simulated  masses.  Exposures  were  made  at 
25kVp,  22mAs  and  a  source-to-image  distance  of  65cm.  Masses  were 
simulated  in  Box  #7,  9,  10, 12.  The  boxes  are  numbered  as  shown  in  Figure  13. 

Film  20  SF  image  of  a  meat  phantom  with  simulated  masses.  Exposures  were  made  at 
25kVp,  22mAs  and  a  source-to-image  distance  of  65cm.  Masses  were 
simulated  in  Box  #7, 9, 10, 12.  The  boxes  are  numbered  as  shown  in  Figure  13. 
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Joint  Photographic  Experts  Group  (JPEG)  Compatible  Data 
Compression  of  Mammograms 

Walter  F.  Good,  Glenn  S.  Maitz,  and  David  Gur 


We  have  developed  a  Joint  Photographic  Experts 
Group  (JPEG)  compatible  image  compression  scheme 
tailored  to  the  compression  of  digitized  mammo- 
graphic  images.  This  includes  a  preprocessing  step 
that  segments  the  tissue  area  from  the  background, 
replaces  the  background  pixels  with  a  constant  value, 
and  applies  a  noise-removal  filter  to  the  tissue  area. 
The  process  was  tested  by  performing  a  just-notice¬ 
able  difference  (JND)  study  to  determine  the  relation¬ 
ship  between  compression  ratio  and  a  reader's  ability 
to  discriminate  between  compressed  and  noncom- 
pressed  versions  of  digitized  mammograms.  We  found 
that  at  compression  ratios  of  15:1  and  below,  image- 
processing  experts  are  unable  to  detect  a  difference, 
whereas  at  ratios  of  60:1  and  above  they  can  identify 
the  compressed  image  nearly  100%  of  the  time.  The 
performance  of  less  specialized  viewers  was  signifi¬ 
cantly  lower  because  these  viewers  seemed  to  have 
difficulty  in  differentiating  between  artifact  and  real 
information  at  the  lower  and  middle  compression 
ratios.  This  preliminary  study  suggests  that  digitized 
mammograms  are  very  amenable  to  compression  by 
techniques  compatible  with  the  JPEG  standard.  How¬ 
ever,  this  study  was  not  designed  to  address  the 
efficacy  of  image  compression  process  for  mammogra¬ 
phy,  but  is  a  necessary  first  step  in  optimizing  the 
compression  in  anticipation  of  more  elaborate  reader 
performance  (ROC)  studies. 

Copyright  ©  1994  by  W.B.  Saunders  Company 

KEYWORDS:  Data  compression,  mammography.  Joint 
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IN  RECENT  YEARS,  there  has  been  a  rapid 
increase  in  the  use  of  mammography  and  a 
corresponding  increase  in  the  complexity  of 
managing  mammographic  studies.  Many  believe 
that  to  accommodate  the  large  number  of  mam¬ 
mographic  studies  now  being  performed,  these 
images  must  be  managed  in  a  digital  format. 
The  main  difficulty  with  this  relates  to  the  large 
volume  of  digital  data  needed  to  represent  each 
image.  High-resolution  digital  mammograms  are 
likely  to  contain  10  to  40  Mbyte  each.  Digital¬ 
archiving  and  data-transmission  systems  quickly 
become  strained  when  called  upon  to  deal  with 
a  large  number  of  these  images.  The  application 
of  image-compression  techniques  has  the  poten¬ 
tial  to  greatly  reduce  this  problem. 

Because  the  problems  of  dealing  with  digital 
image  data  are  universal,  the  Joint  Photo¬ 
graphic  Experts  Group  (JPEG)  has  recently 

Journal  of  DigitaHmaging,  Vol  7,  No  3  (August),  1994:  pp  123-132 


defined  a  general-purpose  image-compression 
standard,^-^  Because  the  JPEG  standard  does 
not  specifically  address  mammographic  images, 
our  aim  in  this  study  was  to  develop  a  compres¬ 
sion  technique  for  digitized  mammograms  which 
is  compatible  with  the  standard,  but  which  is 
optimized  for  mammographic  images,  and  then 
to  determine  the  level  of  compression  that  can 
be  achieved  before  the  effects  of  compression 
become  detectable  by  observers. 

The  JPEG  algorithm  is  a  block  quantization 
version  of  the  cosine  transform  in  which  images 
are  divided  into  8-  x  8-pixel  blocks  and  the 
cosine  transforms  applied  to  each  block  individu¬ 
ally.  The  frequency  coefficients  calculated  from 
the  cosine  transform  are  quantized  by  dividing 
by  values  from  an  array  of  quantization  factors 
and  rounding  the  quotients  to  integral  values. 
The  resulting  array  of  quotients  is  aligned  in  a 
zigzag  order  and  then  encoded  by  either  Huff¬ 
man  coding  or  by  adaptive  arithmetic  coding.  In 
this  method,  the  quantization  factors  that  are 
applied  to  the  frequency  coefficients  determine 
the  compression  ratio  as  well  as  the  kind  of 
information  lost  or  artifact  introduced. 

The  main  challenge  in  applying  JPEG  stan¬ 
dard  compression  relates  to  the  need  to  tailor 
the  technique  to  the  characteristics  of  the  par¬ 
ticular  kind  of  images  being  compressed.  There 
are  two  main  areas  where  there  is  flexibility  in 
applying  the  standard — in  the  quantization  of 
the  frequency  coefficients  and  in  preprocessing 
the  image  before  the  application  of  the  JPEG 
algorithm. 

In  the  past,  most  attention  has  focused  on  the 
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design  of  the  quantization  table  required  in  the 
quantization  stage  of  the  compression  process. 
The  default  quantization  factors,  which  are  pro¬ 
vided  in  the  standard,  are  not  adequate  for  most 
purposes — hence,  the  standard  encourages  user- 
specified  quantization  tables.  But  the  resulting 
image  quality  is  determined  by  the  quantization 
table  used,  so  the  design  of  these  tables  must  be 
optimized  for  the  particular  compression  task. 

The  current  trend  is  to  base  quantization 
tables  on  psychophysical  consideration.^'^  This 
approach  is  somewhat  controversial  when  the 
compression  is  to  be  used  for  radiographic 
images,  because  psychophysical  compression 
techniques  assume  prior  knowledge  of  display 
characteristics  and  viewing  conditions.  The  most 
common  alternative  is  to  use  uniform  quantiza¬ 
tion.  The  justification  for  this  rests  on  the 
assumption  that  there  is  no  a  priori  reason  to 
believe  that  certain  spatial  frequencies  are  more 
important  than  others. 

In  our  interpretation,  the  JPEG  standard 
implies  that  any  compression  process  is  compat¬ 
ible  with  the  standard  if  the  encoded  images  it 
produces  can  be  decoded  with  a  standard  JPEG 
decoder.  Specifically,  this  permits  various  kinds 
of  preprocessing  to  be  applied  to  the  input 
image  data  before  the  actual  compression.  Such 
preprocessing  can  improve  the  overall  results  of 
the  compression/decompression  cycle.  For  ex¬ 
ample,  we  have  shown’  that  when  root  mean 
square  error  (RMSE)  is  used  as  the  measure  of 
image  fidelity,  at  sufficiently  large  compression 
ratios,  there  are  significant  advantages  to  apply¬ 
ing  a  noise-removal  filter  to  the  image  data 
before  compression. 

Another  aspect  of  this  preprocessing,  which  is 
somewhat  unique  to  mammography,  involves 
the  segmentation  of  the  image.  Much  of  the  area 
in  a  typical  digitized  mammogram  corresponds 
to  background  pixels  rather  than  to  tissue  pixels. 
These  background  pixels  contain  noise  that  is 
typically  difficult  to  compress.  To  optimize  com¬ 
pression  in  this  case,  it  is  important  to  minimize 
the  amount  of  data  required  to  represent  these 
background  pixels  in  the  encoded  image.  This 
can  be  done  in  a  JPEG-compatible  manner  by 
first  cropping  the  image  to  the  smallest  rectangu¬ 
lar  area  that  contains  all  of  the  tissue,  and  then 
identifying  the  remaining  background  pixels 
and  replacing  them  with  a  constant  value. 


With  these  considerations  in  mind,  in  our 
laboratory  we  have  developed  and  begun  testing 
software  for  JPEG-compatible  compression  of 
digitized  mammograms. 

MATERIALS  AND  METHODS 

As  described  in  more  detail  below,  our  compression 
software  is  consistent  with  the  extended  12-bit  version  of  the 
JPEG  standard  and  uses  a  set  of  preprocessing  routines 
designed  specifically  for  mammograms.  This  software  was 
written  in  C  to  be  transportable,  but  all  of  the  processing 
reported  in  this  paper  was  performed  on  a  standard  33-MHz 
486  personal  computer.  We  have  tested  this  software  by 
using  it  to  compress  a  set  of  digitized  mammograms  at 
various  levels  of  quantization,  and  then  evaluating  the 
quality  of  these  compressed  images  by  performing  a  just- 
noticeable  difference  (JND)  study. 

Preprocessing.  Before  applying  the  JPEG  compression 
algorithm,  the  image  data  is  preprocessed  to  improve  its 
compressibility.  This  preprocessing  operation  is  divided  into 
two  phases;  first  the  image  is  segmented  and  cropped  and 
the  background  pixel  values  are  modified  to  minimize  the 
storage  requirements  for  the  background,  and  then  a 
noise-removal  filter  is  applied  to  the  tissue  pbcels. 

We  begin  the  segmentation  process  by  automatically 
determining  a  threshold  value  from  the  histogram  of  the 
digitized  source  image  that  indicates  the  transition  between 
background  values  and  tissue  values.  This  threshold  value  is 
used  to  classify  pixels  in  the  image  as  either  tissue  or 
background  so  as  to  produce  a  two-valued  mask  correspond¬ 
ing  to  the  tissue-background  separation.  The  thresholded 
image  normally  consists  of  a  large  contiguous  region  corre¬ 
sponding  to  tissue,  possibly  containing  small  holes,  as  well  as 
smaller  isolated  clusters  of  points  from  the  noisier  back¬ 
ground  regions. 

An  initial  “seed”  point  in  the  tissue  area  is  automatically 
selected  and  grown  to  encompass  the  largest  possible 
connected  region.  This  produces  a  binary  mask  matching 
the  tissue  pbcels  in  the  source  image.  This  mask  is  expanded 
by  15  pixels — ie,  any  pixel  that  is  within  15  pixels  of  a  tissue 
pixel  is  appended  to  the  mask.  This  expansion  is  effected  by 
convolving  the  binary  mask  with  a  31-  x  31-pixel  uniformly 
weighted  kernel  and  then  assigning  all  nonzero  pixels  the 
value  1. 

We  then  convolve  this  binary  mask  with  a  17-  x  17-pixel 
uniformly  weighted  kernel  to  smooth  the  edge  between  the 
mask  (pixels  =  1)  and  background  (pbcels  =  0).  This  ex¬ 
panded  and  smoothed  mask,  whose  pixel  values  represent 
an  interpolation  factor  has  a  value  of  1  for  every 
soft-tissue  pbcel  and  every  background  pbcel  within  7  pixels 
of  a  tissue  pbcel,  and  goes  smoothly  to  zero  over  a  distance  of 
16  pixels  as  the  distance  from  tissue  pbcels  increases.  This 
mask  and  the  source  image  are  both  cropped  to  a  rectangu¬ 
lar  area  that  excludes  as  many  zero  pixels  as  possible  while 
retaining  all  nonzero  pbcels. 

A  constant  value,  which  is  to  be  used  in  modifying 
background  pbcels,  is  determined  by  averaging  the  back¬ 
ground  pbcels  that  are  near  the  skin  boundary.  Finally,  the 
output  image  is  created  pbcel  by  pixel  by  interpolating 
between  the  source  image  and  the  constant  value  using  the 
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formula 

output  =  X  •  source  +  (1  -  X)  •  constant  (1) 

The  resulting  image  has  the  following  characteristics:  (1) 
any  8-  x  8-pixel  block  that  contains  a  tissue  pixel  will  remain 
unchanged  from  the  source  image,  eliminating  the  possibil¬ 
ity  that  the  segmentation  process  will  alter  the  effect  of 
compression  on  any  tissue  pbcel;  and  (2)  background  pixels 
will  have  been  modified  to  change  in  a  smooth  and  continu¬ 
ous  manner  from  their  original  values  to  the  prescribed 
constant  value  as  a  function  of  distance  from  the  nearest 
tissue  pixel,  thus  reducing  the  possibility  that  the  segmenta¬ 
tion  process  will  introduce  an  edge  artifact  at  the  tissue/ 
background  boundary. 

The  second  preprocessing  phase  involves  the  application 
of  a  filter  to  the  tissue  area  of  the  image.  The  particular 
filter  tested  here  is  a  nonlinear  order  filter^*^  that  was  de¬ 
signed  to  modify  mainly  pixels  with  large  noise  components 
while  leaving  the  majority  of  pixels  unchanged.  This  filter 
operates  as  follows.  The  source  image  is  first  decomposed 
into  a  high-frequency  component  and  a  low-frequency 
component  by  unsharp  masking — a  process  whereby  a  low- 
frequency  version  of  the  source  image,  obtained  by  smooth¬ 
ing  the  source,  is  subtracted  from  the  source  image  to 
produce  an  edge-enhanced  (high-frequency)  version  of  the 
source.  Pixel  values  in  each  3-  x  3-pixel  block  of  the 
high-frequency  component  are  sorted.  If  the  center  pbcel  of 
the  block  is  the  highest  value  in  the  block,  it  is  assigned  the 
value  of  the  next  lower  pbcel.  Similarly,  if  it  is  the  lowest 
value,  it  is  assigned  the  next  higher  value.  Otherwise,  the 
pixel  is  left  unchanged.  This  results  in  a  rather  innocuous 
process  of  eliminating  one  pbcel  maxima  or  minima  from  the 
high-frequency  component,  while  leaving  most  pbcels  un¬ 
changed.  The  small  number  of  pixels  changed  are  likely  to 
represent  noise,  as  their  values  are  either  greater  or  less 
than  the  values  of  all  their  immediate  neighbors.  Microcalci¬ 
fications  in  the  image  are  generally  preseved  because  they 
rarely  consist  of  a  single  pbcel.  Finally,  the  filtered  high- 
frequency  component  is  added  back  to  the  low-frequency 
component  to  produce  the  filtered  source  image. 

These  preprocessing  operations  require  two  convolutions 
with  uniform  kernels,  two  applications  of  a  threshold,  a 
region  growing  step,  an  interpolation  step,  and  application 
of  the  nonlinear  filter.  In  our  implementation,  which  does 
not  require  manual  intervention,  all  of  these  steps  are 
performed  in  approximately  one  third  of  the  time  of  the 
compression  itself. 

Image  compression.  The  JPEG  compression  algorithm 
used  in  this  study  is  a  software  package  written  and  tested  in 
our  laboratory  that  implements  the  12-bit  version  of  the 
extended  JPEG  standard.  For  this  study,  Huffman  encoding 
was  used  to  encode  the  quantized  coefficients  and  all 
compression  ratios  were  based  on  this.  The  Huffman  tables 
were  derived  from  the  statistics  of  each  individual  quantized 
image  rather  than  on  the  ensemble  statistics.  Signaling 
information  was  not  included  in  the  code  size  when  calculat¬ 
ing  compression  ratios. 

Results  reported  in  this  paper  are  based  on  uniform 
quantization  of  the  frequency  coefficients.  As  discussed 
below,  any  quantization  scheme  that  .preferentially  de¬ 
grades  the  high-frequency  components  (ie,  most  of  the 


psychophysical  quantization  schemes)  may  be  contraindi¬ 
cated  in  the  compression  of  mammographic  images. 

Evaluation  of  image- compression  technique.  In  an  effort 
to  understand  the  viability  of  the  methods  described  above, 
we  carefully  selected  eight  high-quality  cases  that  spanned 
the  range  of  image  characteristics  customarily  encountered 
in  mammography,  with  some  preference  being  given  to 
images  containing  subtle  pathology.  These  cases  were  digi¬ 
tized  and  the  digital  data  was  compressed  over  a  wide  range 
of  compression  ratios  by  applying  the  above  techniques. 

The  digitization  process  used  a  high-resolution,  high- 
contrast  sensitivity  laser  film  digitizer  (Lumisys,  Sunyvale, 
CA)  that  produces  a  scan  matrbc  of  4,000  x  5,000  pbcels  for 
an  8-  X  10-in  film  by  digitizing  at  a  50-p.m  sampling  interval. 
This  pixel  resolution  results  in  a  Nyquist  spatial  frequency 
of  10  cycles /mm,  which  preserves  the  signal  resolution  of 
the  analog  film.  The  modulation  transfer  function  (MTF)  at 
the  Nyquist  frequency  is  30%  in  the  fast-scan  direction  and 
38%  in  the  slow-scan  direction.  The  16-bit  A'/D  converter  of 
the  digitizer  permits  a  density  measurement  with  an  RMSE 
of  less  than  0.01.  This  represented  less  than  one  third  of  the 
noise  present  in  the  conventional  film  at  densities  exceeding 
0.4. 

Each  digitized  image  was  compressed  at  five  different 
compression  ratios.  Because  it  is  the  degree  of  quantization 
rather  than  the  compression  ratio  itself  that  determines  the 
degradation  of  an  image  by  compression,  the  quantization 
factor  was  used  as  an  independent  variable  and  not  the 
actual  compression  ratio.  The  five  quantization  factors  we 
applied  were  40,  60,  80,  100,  and  120,  which  produced 
average  compression  ratios  between  15:1  and  60:1.  The 
lowest  quantization  factor  is  at  a  level  below  which  it  is  not 
expected  that  one  could  detect  any  effect.  The  highest 
quantization  factor  produced  a  mean  compression  ratio  of 
60:1,  which  is  higher  than  what  is  required  to  permit  the 
efficient  handling  of  mammographic  images  with  current 
digital  technology,  and  it  is  in  the  range  where  studies  of 
other  algorithms  applied  to  other  image  types  have  shown 
significant  deterioration  in  image  quality.’-*^  Note  that  these 
particular  quantization  factors  are  only  appropriate  for 
images  digitized  at  the  resolution  and  with  the  noise 
characteristics  under  consideration  in  this  study. 

All  images  used  in  this  study,  both  the  compressed 
versions  as  well  as  the  original  digitized  data,  were  printed 
on  film  with  a  laser  film  printer  (Eastman-Kodak,  Roches¬ 
ter,  NY).  These  laser-printed  films  were  used  for  all  image 
comparisons.  It  was  not  possible  for  us  to  perform  a 
meaningful  comparison  of  the  laser-printed  images  to  the 
original  (nondigitized)  images  because  our  laser  printer  was 
only  capable  of  printing  at  a  pitch  of  80  fjLm/pbcel,  and  this 
resulted  in  substantial  (1.6  x)  magnification  of  the  laser- 
printed  images.  However,  comparing  two  images  that  have 
both  been  enlarged  to  the  same  degree  is  possible,  and  if  the 
enlargement  had  any  effect,  we  believe  it  would  have  been 
to  make  it  easier  for  readers  to  discriminate  between 
noncompressed  images  and  compressed  images.  The  overall 
digitizer-processor  response  was  calibrated  to  accurately 
duplicate  the  densities  and  contrast  in  the  original  film. 
Noise  contributed  by  the  laser  printer  is  substantially  less 
than  image  noise  caused  by  quantum  mottle  and  film 
granularity  in  the  original  film. 
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Table  1 :  The  Effect  of  Image  Segmentation  on  Tissue  Block 
Representation  (bits  per  block) 


Image 

%  of  Area 
Representing 
Tissue 

Bits  per 
Tissue  Block 
(no 

preprocessing) 

Bits  per 
Tissue  Block 
(after 

segmentation) 

Improvement 

Factor 

1 

51.1 

8.6 

30.1 

3.5 

2 

52.9 

8.5 

27.0 

3.2 

3 

70.5 

9.1 

32.3 

3.5 

4 

71.7 

10.0 

28.0 

2.8 

5 

72.4 

9.4 

23.6 

2.5 

6 

56.8 

7.6 

26.5 

3.5 

7 

70.7 

24.2 

27.4 

1.1 

8 

67.0 

18.0 

22.2 

1.2 

From  the  eight  originals,  which  had  each  been  com¬ 
pressed  at  five  compression  ratios,  a  set  of  40  randomized 
pairs  of  images,  each  consisting  of  one  laser-printed  noncom- 
pressed  image  and  one  compressed  version  of  the  original, 
were  assembled.  Readers  were  asked  to  view  the  pairs  side 
by  side  and  were  forced  to  choose  the  image  they  believed  to 
have  been  compressed.  Readers  were  allowed  to  spend  as 
much  time  as  desired  on  each  pair.  The  reading  environ¬ 
ment  was  equivalent  to  that  normally  used  for  reading 
mammograms  and  readers  were  permitted  to  use  magnify¬ 
ing  lenses. 

Altogether,  20  readers  participated  in  these  readings. 
They  included  radiologists  experienced  in  mammography  as 
well  as  physicists  and  engineers  specializing  in  medical 
imaging. 

RESULTS 

The  performance  of  our  compression  system 
is  somewhat  image  dependent,  but  for  the 
ensemble  of  images  reported  on  here,  we  aver¬ 
aged  about  20  s/Mbyte  to  preprocess  and  com¬ 
press  source  data.  Because  decompression  does 
not  involve  any  steps  similar  to  preprocessing, 
we  were  able  to  perform  it  at  a  rate  of  15 
s/Mbyte  of  output  data.  Although  these  speeds 
are  not  acceptable  for  a  clinical  picture  archiv¬ 
ing  and  communication  system  (PACS)  environ¬ 
ment,  we  believe  they  can  be  easily  improved  by 
adopting  hardware  optimized  for  the  task. 

A  summary  of  the  impact  of  image  segmenta¬ 
tion  on  the  number  of  bits  needed  to  represent 
tissue  blocks  is  presented  in  Table  1.  All  images 
included  in  this  table  had  previously  been 
cropped  to  a  rectangular  area  whose  border  was 
within  1  cm  of  tissue  pixels  in  each  dimension. 
Column  2  gives  the  fraction  of  the  cropped 
image  occupied  by  tissue  pixels.  Columns  3  and 
4  give  the  average  number  of  bits  used  to 
represent  each  8-  x  8-pixel  block  of  tissue  pixels 
before  and  after  application  of  the  segmenta¬ 


tion  process,  for  an  overall  compression  ratio  in 
both  cases  of  40:1.  It  can  be  seen  from  this  table 
that  our  segmentation  process  increased  the 
number  of  bits  used  to  represent  each  tissue 
pixel  by  a  factor  of  2.7  on  the  average  despite 
the  fact  that  the  background  represented  only 
35.9%  of  the  image  area. 

From  approximately  800  observations  in  the 
JND  experiment,  average  observers  correctly 
detected  42%,  55%,  47%,  61%,  and  75%  of  the 
images  for  the  15:1,  23:1,  34:1,  45:1,  and  60:1 
compression  ratios,  respectively.  However,  when 
we  looked  at  the  performance  of  individual 
readers,  we  found  that  there  were  two  readers 
who  were  particularly  good  at  making  the  dis¬ 
criminations  required  in  this  study.  These  read¬ 
ers  were  a  physicist  and  an  electrical  engineer 
who  were  both  specialists  in  image  processing 
and  were  familiar  with  the  particular  artifacts 
caused  by  this  type  of  image  compression.  These 
image-processing  specialists  averaged  discrimi¬ 
nation  rates  of  50%,  67%,  75%,  87%,  and  100% 
for  the  above  compression  ratios.  Their  surpris¬ 
ing  consistency  with  each  other  and  the  mono¬ 
tonic  trend  of  their  performance  led  us  to 
believe  that  we  were  actually  measuring  a  mean¬ 
ingful  threshold.  A  plot  of  the  results  of  this 
JND  study  is  shown  in  Fig  1. 

DISCUSSION 

Performance  of  compression  process.  The  re¬ 
lationship  between  the  compressibility  of  tissue, 
the  compressibility  of  the  background,  and  the 
relative  areas  between  tissue  and  background 


Compression  Ratio 

Fig  1.  Performance  of  image-processing  specialists  and 
nonspecialists  for  the  task  of  identifying  compressed  versus 
noncompressed  images. 
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can  be  expressed  as 

"  A,  +  Ab  _  Ab  ’ 

C  Cb 

in  which  Ct  and  Q  are,  respectively,  the  compres¬ 
sion  ratios  for^^  and^i,,  the  areas  correspond¬ 
ing  to  tissue  and  background,  and  C  is  the 
overall  compression  ratio  for  the  image.  The 
potential  improvement  factor  in  the  number  of 
bits  per  tissue  pixel  that  can  be  achieved  through 
segmentation,  assuming  the  background  is  com¬ 
pressed  to  a  negligible  size  after  segmentation, 
is 

,  l  C 

At  H-  Ab  C5 

where  Q  is  the  compression  ratio  of  the  back¬ 
ground  without  segmentation.  In  all  cases,  after 
segmentation  we  were  able  to  compress  the 
background  at  a  ratio  of  256:1,  which  corre¬ 
sponds  to  using  3  bits  to  represent  each  8-  x 
8-pixel  block.  Theoretically  in  some  unusual 
cases,  it  would  be  possible  to  use  only  2  bits  per 
block,  but  this  is  an  absolute  limit  for  the  JPEG 
algorithm  (sequential  discrete  cosine  transform 
[DCT]  and  Huffman  encoding).  The  actual 
improvements  in  bits  per  tissue  pkel  are  indi¬ 
cated  in  Table  1,  column  5. 

Note  that  images  7  and  8  behaved  somewhat 
differently  than  the  other  six  images.  The  back¬ 
grounds  of  these  images  were  so  dense  that  they 
were  beyond  the  range  of  our  digitizer.  Conse¬ 
quently,  many  background  pixels  values  were 
saturated  at  the  maximum  output  of  our  digi¬ 
tizer,  and  hence,  the  apparent  noise  in  the 
background  was  greatly  reduced.  Saturation  of 
the  background  had  much  the  same  effect  as  our 
segmentation  process,  but  has  the  potential  to 
also  alter  tissue  pixels  in  an  undesirable  way. 
We  have  included  these  particular  images  after 
verifying  that  the  tissue  pixels  themselves  were 
not  saturated. 

For  images  1  through  6,  our  segmentation 
process  increased  the  number  of  bits  used  to 
represent  each  tissue  pixel  by  a  factor  of  3.2  on 
the  average.  If  images  7  and  8  are  also  included, 
this  effect  is  reduced  to  2.7.  The  average  area 
occupied  by  tissue  pfacels  in  cropped  images  1 


through  6  was  62.6%.  It  can  be  seen  that  for 
images  in  which  the  background  is  not  satu¬ 
rated,  the  average  improvement  in  bits  per 
tissue  pixel  is  much  larger  than  what  would  be 
suggested  simply  by  the  relative  areas  occupied 
by  the  tissue  pixels  and  background  pixels. 

We  illustrate  what  happens  during  the  seg¬ 
mentation  process  by  using  image  number  6 
from  our  study  as  a  specific  example.  This 
image,  which  is  partially  shown  in  Fig  2,  is 
typical  of  what  we  have  observed  in  general. 
Initially,  this  mammogram  was  digitized  to  a 
12.3-Mbyte  file.  Cropping  reduced  this  by  5,7% 
to  11.6  Mbyte,  and  of  this  file,  56.8%  of  the 
pixels  corresponded  to  tissue. 

A  plot  of  the  compressibility -of  various  re¬ 
gions  of  this  image,  which  are  identified  in  Fig  2, 
is  presented  in  Fig  3.  The  independent  variable 
on  this  graph,  the  quantization  factor,  is  speci¬ 
fied  at  the  time  of  compression  and  directly 
controls  the  loss  of  information  during  the 
compression  process.  It  can  be  seen  that  differ¬ 
ent  regions  of  the  image  vary  widely  with  re¬ 
spect  to  their  compressibility  and  that  the  back¬ 
ground  is  significantly  less  compressible  than 
either  of  the  tissue  regions.  In  many  cases  such 
as  this  one,  because  the  background  is  very 
noisy,  it  is  the  least  compressible  part  of  the 
image.  When  the  nonsegmented  version  of  this 
image  was  compressed  at  a  ratio  of  25:1  (0.48 
bits  per  pixel),  the  compression  ratio  of  the 
background  was  only  14.9:1  compared  with 
51.8:1  (0.23  bits  per  pixel)  for  the  tissue  pixel. 
The  encoded  file  contained  0.464  Mbyte  of  data 
of  which  0.337  Mbyte  was  needed  to  describe 
the  background  and  the  remaining  0.127  Mbyte 
represented  tissue.  In  other  words,  despite  the 
fact  that  the  background  was  only  43%  of  the 
source  image,  more  than  72%  of  the  data  in  the 
encoded  file  was  used  to  describe  the  back¬ 
ground.  After  the  image  was  preprocessed  and 
compressed  at  the  same  25:1  compression  ratio, 
the  constant  background  in  the  segmented  im¬ 
age  was  now  compressed  at  a  ratio  of  256:1  (.05 
bits  per  pixel)  and  the  tissue  area  was  com¬ 
pressed  at  a  ratio  of  14.8:1  (0.81  bits  per  pixel). 
The  background  was  reduced  to  only  0.02  Mbyte 
or  about  4%  of  the  encoded  data.  The  net  effect 
was  that  segmentation  increased  the  amount  of 
data  used  to  represent  tissue  pixels  by  a  factor 
of  3.5.  Figure  4  compares  the  overall  compress- 
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Fig  2.  image  showing  the 
three  typical  regions  of  interest 
in  different  areas  representing 
tissue  and  background. 


ibility  of  the  segmented  image  to  that  of  the 
nonsegmented  image.  The  compressibility  of  a 
chest  image  is  also  included  for  comparison. 

Figure  5  shows  in  more  detail  the  relative 
number  of  bits  per  tissue  pixel  between  the 


Fig  3.  Comparison  of  the  compressibiilty  of  the  three 
regions  marked  in  Fig  2,  with  the  compressibility  of  the  image 
as  a  whole. 


segmented  and  nonsegmented  images.  The  ac¬ 
tual  pattern  of  the  compressibility  of  various 
parts  of  the  image  is  also  of  interest  because  it 
provides  information  that  eventually  will  be 
needed  for  optimizing  the  relative  compression 


Fig  4.  Comparison  of  the  compressibility  of  mammograms 
with  and  without  preprocessing  to  the  compressibility  of 
digitized  chest  images. 
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Fig  5.  Comparison  of  the  mean  number  of  bits  used  to 
represent  a  tissue  pixel  for  segmented  versus  nonsegmented 
images. 


of  different  tissues.  Figure  6,  A  and  B  show  the 
pattern  of  bits  per  pixel  with  and  without 
preprocessing.  It  is  clear  from  a  comparison  of 
these  images  that  preprocessing  results  in  a 
dramatic  shift,  from  bits  being  used  to  represent 
background  to  being  used  to  represent  tissue. 
One  further  observation  from  Fig  6B  is  that  a 
large  fraction  of  the  bits  used  to  represent  tissue 
are  in  fact  used  to  represent  only  the  regions 
near  skin.  This  observation  is  consistent  through¬ 
out  the  ensemble  of  images  we  have  studied  to 
date.  We  believe  it  is  mostly  caused  by  the  low 
compressibility  of  the  high-frequency  structures 
in  the  near-skin  regions  and  to  a  lesser  extent 
caused  by  the  high  density  of  the  skin  areas  on 
the  original  films  that  produces  increased  digi¬ 
tizer  noise.  This  can  be  seen  in  Fig  7,  which 


Fig  6.  Images  of  bits  per  block 
before  segmentation  (A)  and  af¬ 
ter  segmentation  (B)  for  an  im¬ 
age  compressed  at  40:1. 
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relates  the  number  of  bits  per  pixel  to  film 
density  for  a  fixed  quantization  level,  for  reason¬ 
ably  flat  areas  of  an  image.  The  increasing  trend 
as  a  function  of  optical  density  is  caused  by  the 
increasing  noise  of  the  digitizer  as  film  density 
increases.  Also  shown  in  Fig  7  are  measure¬ 
ments  of  the  compressibility  of  skin  areas  as  a 
function  of  density.  These  points  are  clearly 
above  the  line.  We  have  observed  that  this  low 
compressibility  of  the  skin  areas  also  holds  for 
storage  phosphor  (computed  radiography)  im¬ 
ages  despite  the  very  different  noise  characteris¬ 
tics  of  these  images. 

Future  updates  to  the  JPEG  standard  are 
likely  to  contain  provisions  for  adaptive  quanti¬ 
zation.  When  this  happens,  it  will  be  necessary 
to  address  the  issue  of  the  relative  number  of 
data  bits  that  should  be  allocated  to  represent 
each  particiular  type  of  tissue. 

Ejfect  of  filtration.  The  effects  of  applying 
the  nonlinear  filter  to  mammogram  number  6  is 
shown  in  Fig  8.  Although  the  filter  only  changed 
6.7%  of  the  pixel  values,  it  seemed  to  provide  a 
net  benefit  at  the  higher  compression  ratios.  Its 
behavior  was  comparable  with  what  was  achieved 
when  a  similar  filter  was  applied  to  the  compres¬ 
sion  of  chest  images.’  Figure  8  shows  the  RMSE 
for  the  compression  of  the  filtered  source  along 
with  the  RMSE  for  the  compression  of  the 
nonfiltered  source  as  they  depend  on  the  com¬ 
pression  ratio.  It  can  be  seen  that  at  low 
compression  ratios,  the  RMSE  of  the  filtered 
image  is  higher  because  of  changes  in  the  image 
caused  by  the  filter,  but  at  higher  compression 


Fig  7.  Plot  of  compressibility  of  different  ''nonskin"  tissue 
regions  as  a  function  of  film  density  (line).  Also  included  for 
comparison  are  data  from  different  points  in  the  skin  region. 


Fig  8.  RMSE  resulting  from  nonlinear  filtration  as  a  function 
of  compression  ratio.  Results  for  both  the  filtered  and  nonfil¬ 
tered  image  are  shown. 

ratios,  the  errors  caused  by  compression  over¬ 
whelm  the  error  caused  by  filtration  and  be¬ 
come  predominant.  The  compression-induced 
errors  are  consistently  smaller  for  the  filtered 
image  than  for  the  nonfiltered  image,  and  at 
higher  compression  ratios,  the  filtered  image 
actually  deviates  less  from  the  source  image. 
This  behavior  persisted  across  our  image  set, 
but  was  found  to  be  image  dependent.  It  is  our 
view  at  this  time  that  filtration  of  this  sort  is 
advantageous  at  compression  ratios  above  20:1, 
but  it  remains  to  be  tested  in  diagnostic  perfor¬ 
mance  studies. 

Ejfect  of  segmentation  and  cropping.  In  gen¬ 
eral,  for  all  of  our  images,  preprocessing  the 
digitized  data  made  a  dramatic  improvement  in 
compressibility.  This  improvement  can  be  attrib¬ 
uted,  for  the  most  part,  to  the  benefit  of  segmen¬ 
tation.  It  is  actually  the  degree  of  compression 
of  the  image  area  corresponding  to  tissue  pixels, 
or  the  number  of  bits  per  tissue  pixel,  rather 
than  the  compression  ratio  for  the  entire  image 
array,  that  directly  determines  the  impact  of 
compression  on  the  diagnostic  quality  of  mam- 
mographic  images.  In  mammography,  there  is 
wide  variation  from  case  to  case  in  the  propor¬ 
tion  of  the  film  occupied  by  tissue  pixels.  In  all 
cases,  after  segmentation  and  cropping,  the 
background  will  be  compressed  at  a  very  high 
compression  ratio,  whereas  the  tissue  pixels  are 
compressed  at  a  more  modest  ratio. 

If  images  are  not  segmented,  then  cropping  is 
particularly  beneficial.  However,  if  the  images 
are  segmented  by  the  above  procedure,  the 
benefit  of  cropping  is  largely  neutralized  by  the 
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efficiency  of  the  segmentation  process  in  increase 
ing  the  compressibility  of  the  background,  Nev- 
ertheiess,  cropping  can  be  useful,  because  when 
the  image  is  reconstructed,  it  will  retain  its  size 
before  encoding.  Maintaining  excessive  back¬ 
ground  pixels  could  have  adverse  effects  on  the 
time  required  to  encode  and  decode  the  image 
as  well  as  on  other  system  components  such  as 
image  displays. 

Preservation  of  image  quality.  Much  of  a 
mammographic  image  does  not  contain  appre¬ 
ciable  high-frequency  information,  but  the  need 
to  preserve  microcalcifications  with  dimensions 
of  less  than  100  jim,  which  sometimes  appear  in 
mammographic  films,  may  require  that  we  digi¬ 
tize  mammograms  at  50  pm.  This  has  the  effect 
of  oversampling  most  of  the  image — ie,  the 
actual  information  content  of  the  digitized  mam¬ 
mogram  is  much  less  than  the  size  of  the  image 
matrix  would  suggest. 

It  does  not  require  severe  quantization  to 
achieve  a  high  compression  ratio  on  areas  that 
do  not  contain  high-frequency  information.  The 
remainder  of  the  image,  such  as  areas  contain¬ 
ing  microcalcifications  or  the  skin,  is  not  as 
compressible,  but  these  areas  generally  com¬ 
prise  such  a  small  fraction  of  the  total  area  that 
they  have  minimal  impact  on  the  compression 
ratio.  Together,  these  factors  imply  that  it  may 
not  be  necessary  or  desirable  to  preferentially 
degrade  the  higher-frequency  components  in 
the  quantization  process  to  achieve  high  com¬ 
pression  ratios,  and  consequently,  it  is  possible 
to  preserve  features  such  as  microcalcifications 
during  the  compression  process. 

In  comparing  the  compressibility  of  mammo¬ 
grams  to  that  of  chest  images,  we  found  that 
preprocessed  mammograms  are  very  similar  to 
chest  images  in  compressibility;  however,  un¬ 
processed  mammograms  are  much  more  diffi¬ 
cult  to  compress  than  chest  images.  We  have 
included  a  plot  in  Fig  4  that  was  derived  from  a 
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previous  study’  of  the  compression  ratio  of 
chest  images  versus  the  quantization  factor.  The 
chest  image  had  been  digitized  at  100-p.m  reso¬ 
lution  with  12  bits  of  dynamic  range.  At  each 
level  of  quantization,  the  RMSE  for  the  chest 
image  was  similar  to  the  RMSE  for  mammo¬ 
grams  compressed  with  the  same  quantization 
factor. 

Observer  performance.  When  the  differences 
between  images  are  small,  it  is  possible  that 
readers  can  detect  the  difference,  but  still  not 
be  able  to  tell  which  image  is  compressed.  The 
key  to  differentiating  between  images  at  the 
lower  compression  ratios  seemed  to  be  the 
presence  or  absence  of  blocking  artifact.  For  the 
compression  techniques  we  applied  in  this  study, 
blocking  artifact  becomes  noticeable  before  any 
degradation  of  image  resolution  is  obvious. 
Blocking  artifact  could  be  detected  by  image- 
processing  specialists  and  nonspecialists  alike, 
but  the  nonspecialist  readers  were  more  likely 
to  misinterpret  it.  It  is  clear  from  discussions 
with  some  readers  upon  completion  of  their 
readings,  that  when  the  blocking  artifact  was 
barely  visible  they  sometimes  mistook  it  for 
actual  image  information  and  assumed  the  im¬ 
age  with  the  artifact  was  the  noncompressed 
image.  The  two  best  readers  seemed  to  be  able 
to  avoid  this  mistake. 

The  JND  study  presented  here  is  not  in¬ 
tended  to  address  the  question  of  the  impact  of 
compression  on  diagnostic  performance.  Its  pur¬ 
pose  was  merely  to  give  us  an  opportunity  to  test 
the  feasibility  of  the  techniques  we  believe  are 
appropriate  for  the  compression  of  mammo¬ 
graphic  data,  and  to  identify  the  range  of 
compression  ratios  that  merit  further  study. 
There  are  many  possibilities  for  further  refine¬ 
ments,  particularly  in  the  area  of  preprocessing 
the  data,  that  are  within  the  constraints  of  the 
standard,  and  these  will  likely  extend  the  useful 
range  of  compression  ratios. 
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